UNCLA 

24  365  AN  EXPERIMENTAL  STUDY  OP 
TURBULENT  BOUNDARY  LAVER 
THERHOSCIENCES  DIV  R  JA 
SSIFIED  ARO-16760  2-EG  ARO-1-80 

THE  DYNAMICS 
<U>  STANFORD 
VARAHAN  ET  AL 

OF  AK 
UNIV  C 
DEC 

UNSTEADY 

A 

82  TF-18 
.  F/G  28/4 

n 

■ 

"H5SB 

HHfll 

a 

HIH 

ADA  124365 


AN  EXPERIMENTAL  STUDY  OF  THE  DYNAMICS 
OF  AN  UNSTEADY  TURBULENT  BOUNDARY  LAYER 


by 

Rangarajan  Jayaraman 
Pradip  Parikh 
William  C.  Reynolds 


Prepared  under  the  joint  support  of 
Army  Research  Office 

Aeromechanics  Laboratory,  RTL  (AVRADCOM) 
NASA/Ames  Research  Center 


Technical  Report  Number  TF-18 


Thermosciences  Division 
Department  of  Mechanical  Engineering 
Stanford  University 
Stanford,  California 


December  1982 

83  02  014  07  4 


Unci  ass  i  f  ied _ /'T'T  f  // y 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fWTi«n  Doto  Snltrtd)  \^  /  ! 

REPORT  DOCUMENTATION  PAGE  b eforeP compl  etin  g*form 

1.  REPORT  NUMBER  |2.  GOVT  ACCESSION  NO.  3-  RECIPIENT’S  CATALOG  NUMBER 


REPORT  DOCUMENTATION  PAGE 


1.  REPORT  HUMBER 

16760. 2-EG 

4.  TITLE  (and  Submit) 


An  Experimental  Study  of  the  Dynamics  of  an 
Unsteady  Turbulent  Boundary  Layer 

7.  AUTHOR/..) 

Rangarajan  Jayaraman 

Pradip  Parikh 

V/i  1 1  iam  C.  Reynolds 

V.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Stanford  University 
Stanford,  CA 


5.  TYPE  OF  REPORT  4  PERIOD  COVERED 

Final : 

1  Oct  79  -  30  Sep  82 

6.  PERFORMING  ORG.  REPORT  NUMBER 
4.  CONTRACT  OR  GRANT  HUMBERT*) 

ARO  1-80;  2-81;  1-32 


10.  PROGRAM  ELEMENT.  PROJECT,  TASK 
AREA  4  WORK  UNIT  NUMBERS 


II.  CONTROLLING  OFFICE  NAME  AND  AOORESS 

U.  S.  Army  Research  Office 

Post  Office  Box  12211 

Research  Triangle  Park,  NC  27709 

14.  MONITORING  AGENCY  NAME  4  ADDRESS///  dll  It  ton  I  i 


I*.  REPORT  DATE 


l>.  NUMBER  OF  PAGES 


Controlling  Offtct)  I  IS.  SECURITY  CLASS,  (ol  Uilo  raporfj 


€ 


16.  DISTRIBUTION  STATEMENT  (ol  thlo  Roport) 


IS*.  DECLASSIFICATION/ DOWNGRADING 
SCHEDULE 


Approved  for  public  release;  distribution  unlimited. 


r\"HO 

9  0  ^*6*4# 


17.  DISTRIBUTION  STATEMENT  (of  tho  •6«fr«cf  mltrtd  In  Bloc*  JO,  If  different  Ami  ft«pd| ^ 


V 


T51? 


!  It*  supplementary  notes 


The  View,  opinions,  and/or  findings  contained  in  this  report  are  those  of  the 
author (s)  and  should  not  be  construed  as  an  official  Department  of  the  Army 
position,  policy,  or  decision,  unless  so  designated  by  other  documentation 

19.  KEY  WORDS  (Continue  on  r#v«rM  tide  If  nocoemory  «id  identify  by  block  number) 

turbulent  boundary  layers 

fluid  dynamics 

free  stream  velocity 

B4  ABSTRACT  fCmattmm  an  roro too  oSS  II  rmoomoory  and  Idontitr  4y  block  mi6*rj 

VlTurbulent  boundary  layers  with  a  time-dependent  free-stream  velocity  occur  in 
many  fluid  dynamic  devices.  The  main  objectives  of  the  present  work  are  as 
follows:  to  understand  the  characteristics  of  a  turbulent  boundary  layer  that 
develops  under  steady,  zero-pressure  gradient  conditions  upstream,  and  Is  sub¬ 
jected  to  well-defined  and  sinusoidally  varying  adverse  free-stream  gradient 
downstream,  and  to  provide  a  definitive  data  set  that  would  serve  as  a  test  cas< 
for  evaluating  the  performance  of  calculation  procedures  and  turbulence  models  or 
such  flows.*—' 


FORM 

i  jam  n 


tTIOM  OF  I  MOV  66  16  OBSOLETE 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fW*r*n  Doto  Bntotod) 


AN  EXPERIMENTAL  STUDY  OF  THE  DYNAMICS 
OF  AN  UNSTEADY  TURBULENT  BOUNDARY  LAYER 


by 

Rangarajan  Jayaraman 
Pradip  Parikh 
William  C.  Reynolds 


Prepared  under  the  joint  support  of 


Army  Research  Office 
Aeromechanics  Laboratory,  RTL  (AVRADCOM) 
NASA/Ames  Research  Center 


Technical  Report  No.  TF-18 


Thermosciences  Division 
Department  of  Mechanical  Engineering 
Stanford  University 
Stanford,  California  94305 


December  1982 


Ac knowl edgments 

We  gratefully  acknowledge  the  continued  sponsorship  of  the  Army 
Aeromechanics  Laboratory,  NASA-Ames  Research  Center,  and  the  Army 
Research  Office.  Dr.  Larry  Carr  of  the  Amy  Research  Office  attended 
all  our  weekly  meetings  and  provided  us  with  many  valuable  suggestions. 
We  thank  Dr.  Carr  and  Prof.  Stephen  Kline  for  carefully  reviewing  the 
manuscript  and  suggesting  many  improvements. 

Messrs.  Robin  Birch,  Tots  Ikebe,  and  Roy  Niegel  provided  excellent 
shop  support  throughout  this  work.  Mrs.  Ruth  Korb  has  done  an  outstand¬ 
ing  job  with  all  phases  of  the  preparation  of  this  report.  Her  efforts 
and  editorial  expertise  are  very  much  appreciated. 


jrlr 


~iT»u«uu. 


Abstract 


Turbulent  boundary  layers  with  a  time- dependent  free- stream  veloc¬ 
ity  occur  In  many  fluid  dynamic  devices.  The  main  objectives  of  the 
present  work  are  as  follows: 

•  to  understand  the  characteristics  of  a  turbulent  boundary 
layer  that  develops  under  steady,  zero-pressure  gradient 
conditions  upstream,  and  is  subjected  to  well-defined  and 
sinusoidally  varying  adverse  free- stream  gradient  down¬ 
stream,  and 

•  to  provide  a  definitive  data  set  that  would  serve  as  a 
test  case  for  evaluating  the  performance  of  calculation 
procedures  and  turbulence  models  for  such  flows. 

The  experiments  were  done  in  a  unique  closed-loop  water  channel 
especially  constructed  for  the  study.  A  single-color,  forward- scatter, 
frequency- shifted,  Laser-Ooppler  Anemometer  system  was  used  with  a 
frequency  tracker  to  measure  the  streamwise  component  of  velocity.  A 
laboratory  microcomputer  system  was  set  up  to  control  the  experiment  and 
automate  the  data-acqulsltlon  process.  Two  major  series  of  experiments 
with  relative  amplitudes  (of  free-stream  oscillations  at  the  end  of  the 
test  section)  of  about  5%  and  25%  were  done.  Several  frequencies  and 
X  stations  were  covered  in  each,  to  provide  a  wide  range  in  Strouhal 
number  (0.04  <  S6  <  1.6).  Phase- averaged  streamwise  velocity  and 
square  of  streamwise  turbulent  fluctuations  were  measured;  from  these, 
other  quantities  of  interest  were  computed. 

The  results  show  that  the  time-averaged  velocity  and  turbulence  are 
insensitive  to  the  imposed  unsteadiness,  while  the  periodic  components 
are  very  strongly  affected  by  the  frequency.  The  periodic  velocity 
field  exhibits  two  distinct  regimes  of  behavior  with  respect  to  the 
Strouhal  number  (Sx,  based  on  the  local  mean  free-stream  velocity  and 
the  streamwise  distance  in  the  test  section).  For  Sx  <  1,  the  flow  is 
primarily  governed  by  Inertia,  while  for  Sx  >  1,  the  flow  Is  dominated 
by  Increasing  pressure-gradient  oscillations  that  overcome  inertia.  The 
high  Strouhal  number  results  agree  quite  well  with  the  Stokes  solution. 
The  steady  form  of  the  law  of  the  wall  seems  to  be  Inapplicable  to  the 
phase- averaged  velocity  profiles  at  high  frequencies.  Uver  the  range  of 
frequencies  and  amplitudes  studied,  there  seem  to  be  no  dramatic  Inter¬ 
actions  between  the  organized  oscillations  and  the  disorganized  turbu¬ 
lent  fluctuations. 
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Nomenclature 

A  First  harmonic  amplitude  of  free-stream  velocity  oscillations; 

used  In  order-of-magnltude  analysis. 

A^f  nth  harmonic  amplitude  of  variable  ft  unless  noted  otherwise. 

Al,01  First  harmonic  amplitude  of  <6j>;  used  In  plots. 

A1,U  *  A1>u;  used  Pl°ts* 

ai,ue)  .  .  .  , 

j  *  Al  u  J  used  1n  Plots. 

Al,Ue  J  *  * 

A1,U'U‘  *  Al,u'u'i  used  1n  Plots 

Al,uT  F^st  harmonic  amplitude  of  <  ut  >• 

a  s  A1,U  /uo  evaluated  at  X*  «  1;  first  harmonic  amplitude  of 

free-stream  velocity  at  the  end  of  the  test  section  normalized 
by  the  mean  free-stream  velocity  at  the  entrance  to  the  test 
section. 

C  *  5.0;  constant  In  log  law  of  the  wall. 

CF  *  Cf;  used  In  plots. 

Cp  «  Cp;  used  In  plots. 

Cf  *  Zrw/pu|  ;  skin-friction  coefficient  of  the  mean  velocity 

profile  determined  by  Coles  analysis  procedure;  not  an 
Instantaneous  quantity. 

2 

Cp  *  2(p-p0)/pU‘;  coefficient  of  pressure. 

2 

Cp  *  2tPs‘po^f*Jo;  stead*  stattc  pressure  coefficient, 
s 

Dc  »  6;  used  In  plots. 

D  *  Uc;  not  a  time  average. 

Dgg  Mean  value  of  <6gg>;  used  In  plots. 

0}  *  6j;  used  In  plots. 

Dj  Mean  value  of  <  6}  >;  used  In  plots. 

Dg  ■  used  In  plots. 

<0g9>  -  <2gg>;  used  In  plots. 
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■  <6}>;  used  In  plots. 

*  <$2>;  used  In  plots. 

Frequency  In  Hz;  also  any  flow  variable 


<DX> 

<02> 

f 

S  »  dCp/dX';  normalized  streamwlse  gradient  of  coefficient  of 

pressure. 

4  Shape  factor  of  the  mean  velocity  profile;  not  an  Instantaneous 

quantity. 

CH>  Shape  factor  computed  from  phase-averaged  velocity  profiles; 

not  a  phase- averaged  quantity;  used  In  plots. 

1  Mean  value  of  <H>;  used  In  plots. 

i  Index  In  subscript  notation. 

j  Pure  Imaginary  number;  also  used  as  a  second  Index  In  subscript 

notation. 

L  Length  of  the  test  section,  61  cm;  also  typical  streamwlse 

length  scale  used  In  order-of-magnltude  analysis. 

L«j  «  J^;  used  In  plots. 

PL  Left-hand  side  of  normalized  momentum  Integral  equation. 

Pj_  *  ^1,U  at  the  lest  measured  X  station  In  recovery  section; 

used  In  plots. 


pl, dCp/dX'  *  *1 .dCp/dX1 1  used  1n  P1ots* 

P1,D1  Flrst  harmonic  phase  of  <  6j  >;  used  In  plots. 

1,U  *  *l,u;  used  1°  plots. 

1,UEJ  ■  to  n  ;  used  In  plots. 
l,Ue)  A,u« 


1,U'U' 


*  ^i.u'u1*  used  In  plots. 

Right-hand  side  of  normalized  momentum  Integral  equation. 
Instantaneous  static  pressure. 

Reference  static  pressure  at  the  entrance  to  the  test  section. 
Steady  static  pressure 


*  U  6g/v;  momentum  thickness  Reynolds  number  for  mean  velocity 
profile;  used  In  plots. 

*  uL/U  ;  Strouhal  number  used  In  order-of-magnltude  analysis 

e 

based  on  the  typical  length  scale  In  the  streamwlse  direction 
and  the  mean  free- stream  velocity. 

*  <o6/Ue  ;  Strouhal  number  based  on  local  values  of  the  thick¬ 
ness  of  the  mean  boundary  layer  and  the  free-stream  mean  veloc¬ 
ity. 

8  co(X-X  )/U  ;  Strouhal  number  based  on  the  streamwlse  distance 

^  o'  e 

In  the  test  section  at  the  measuring  station  and  the  local  mean 
free-stream  velocity. 

Time. 

Discrete  times  at  which  a  variable  Is  sampled. 

■  U;  time  average  of  streamwlse  velocity. 

8  U  ;  used  In  order-of-magnltude  analysis, 
e 

Uniform  mean  free-stream  velocity  In  the  upstream  development 
section;  approximately  730  mm/s. 

8  U/u  . 
x 

8  u'u';  used  In  plots. 

8  uT;  used  In  plots. 

8  u;  used  In  plots. 


8  <u>. 

Instantaneous  streamwlse  component  of  velocity. 

1  8  1,2,3;  components  of  the  velcolty  vector  In  subscript 
notation. 

Wall  shear  velocity  of  the  mean  velocity  profile  computed  by 
Coles  analysis  procedure;  not  an  Instantaneous  quantity. 

"Wall  shear  velocity"  parameter  used  In  the  log-linear  region 
analysis;  not  necessarily  related  to  the  phase-averaged  wall 
shear  stress. 

8  v;  time  average  of  normal -to-wall  velocity. 

Instantaneous  normal -to-wall  component  of  velocity. 
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SB 


M  Channel  width,  30.48  cm. 

w  Spanwlse  component  of  velcolty. 

X  St reamwlse  distance  In  the  experimental  rig  measured  with  re¬ 

spect  to  the  nozzle  exit. 

XQ  Streamwlse  distance  of  the  entrance  to  the  test  section  with 
respect  to  the  nozzle  exit,  205.1  cm. 

X'  »  (X-X0)/L;  streamwlse  distance  of  the  measuring  station  with 

respect  to  the  entrance  to  the  test  section  normalized  by  the 

test-section  length. 

x  *  X. 

jc  Position  vector  In  30  space. 

Xj  1  -  1,2,3;  components  of  the  position  vector  In  subscript 

notation. 

Y  l^nnal-to-wall  distance  In  the  experimental  rig  measured  with 

respect  to  the  wall. 

Y+  -  yu^/v. 

ys  *  y/Jtjj  normal -to-wall  distance  normalized  by  Stokes  layer 

thickness  parameter. 

y  -  Y. 

Z  Spanwlse  distance  In  the  experimental  rig  measured  with  respect 

to  the  centerline  of  the  channel. 

z  -  Z. 


a  •  A/U e;  first  harmonic  amplitude  of  free- stream  velocity 

normalized  by  the  mean  free-stream  velocity;  used  In  order-of- 
magnltude  analysis. 


6  Thickness  of  the  mean  velocity  profile;  usually  determined  by 

Coles  analysis  procedure  (see  Coles  [1968]). 

6fN  N^-order  uncertainty  In  variable  f. 


<6gg>  99%  thickness  computed  from  phase- averaged  velocity  profiles; 

not  a  phase-averaged  quantity. 
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Displacement  thickness  of  the  mean  velocity  profile;  not  an 
Instantaneous  quantity. 

Displacement  thickness  computed  from  phase-averaged  velocity 
profiles. 

Typical  length  scale  In  the  normal -to-wall  direction  for  peri¬ 
odic  quantities  used  In  the  order-of-magnltude  analysis;  not  a 
turbulent  quantity. 

Momentum  thickness  of  the  mean  velocity  profile;  not  an  Instan¬ 
taneous  quantity. 

Momentum  thickness  computed  from  phase-averaged  velocity 
profiles;  not  a  phase- averaged  quantity. 

Multiplicative  factor  used  In  order-of-magnltude  analysis. 


Von-Karman  constant  used  In  log  law  of  the  wall.  U.41. 

r STv 

■If —  ;  Stokes  layer  thickness  parameter, 
v  w 

Kinematic  viscosity. 

Circular  frequency  In  radians/ sec. 

Phase  of  nth  harmonic  component  of  periodic  variable  f. 
First  harmonic  phase  of  <u,>. 

Density. 


*  2n/u>;  time  period  of  Imposed  excitation. 
Wall  shear  stress. 

Cycle  phase  angle. 


il  Symbols  and  Usage 

nth  harmonic  amplitude  of  variable  f  unless  noted  otherwise. 

Phase  of  n^  harmonic  component  of  periodic  variable  f  unless 
noted  otherwise. 


Used  over  a  variable  denotes  Its  time  average,  unless  noted 
otherwise. 


Used  over  a  variable  denotes  Its  periodic  component. 

Used  after  a  variable  denotes  Its  typical  value  of  order-of- 
magnltude  accuracy. 

Used  as  a  superscript  to  a  variable  denotes  Its  turbulent 
component  unless  noted  otherwise. 

Used  as  a  subscript  to  a  variable  denotes  Its  free-stream 
val ue. 

Used  under  a  variable  denotes  that  It  Is  a  vector. 

Used  before  a  variable  denotes  changes  In  that  variable. 

Used  around  a  variable  denotes  Its  phase  average  unless  noted 
otherwise. 

Used  as  a  subscript  to  a  variable  denotes  Its  root-mean- squared 
value. 
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Chapter  1 
INTRODUCTION 


1.1  General 

Turbulent  boundary  layers  with  a  time-dependent  (unsteady)  free- 
stream  velocity  occur  In  many  fluid  dynamic  devices.  Examples  of  such 
flows  are  flow  over  the  blades  of  a  helicopter  In  forward  flight,  flow 
over  the  wings  of  an  accelerating  aircraft,  and  flow  over  the  stator 
blades  downstream  of  a  rotary  stage  In  a  turbomachine. 

The  work  described  in  this  report  deals  with  an  experimental  study 
of  the  dynamics  of  a  turbulent  boundary  layer  subjected  to  an  oscilla¬ 
tory,  mean  adverse,  free-stream  velocity  gradient.  It  Is  part  of  an 
ongoing  research  program  at  Stanford  Investigating  the  effects  on  tur¬ 
bulence  of  Imposed  organized  unsteadiness. 

There  are  two  major  objectives  for  the  present  study.  First,  the 
work  seeks  to  provide  an  understanding  of  the  physics  of  unsteady  turbu¬ 
lent  boundary  layers.  Second,  the  experiments  are  intended  to  provide  a 
definitive  data  base  for  the  evaluation  and  development  of  turbulence 
models  for  such  flows. 

To  provide  further  motivation  and  background  for  the  present  work, 
a  brief  review  of  previous  research  which  Is  believed  to  be  most  rele¬ 
vant  Is  given  In  the  next  section. 

1.2  Previous  Work 

Table  1.1  Is  a  summary  of  Important  parameters  and  pertinent  infor¬ 
mation  from  several  of  the  experiments  reviewed  below.  Experiments  on 
unsteady,  fully  developed  channel  and  pipe  flows  are  also  Included, 
since  the  near-wall  turbulence  structure  In  such  flows  Is  similar  to  the 
turbulent  boundary  layer. 

Two  Important  parameters  Included  In  the  table  are  an  amplitude 
parameter  which  Is  the  free-stream  oscillation  amplitude  normalized  by 
the  mean  free-stream  velocity  (a)  and  a  frequency  parameter  In  the 
form  of  a  Strouhal  number  based  on  the  boundary  layer  thickness  and  the 
mean  free-stream  velocity  (S^).  For  fully  developed  pipe  flows,  the 
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rad t us,  and,  for  fully  developed  channel  flows,  the  channel  half-gap  are 
used  In  place  of  the  boundary  layer  thickness;  the  mean  centerline  vel¬ 
ocity  Is  used  In  place  of  the  mean  free-stream  velocity.  Some  of  the 
parameters  had  to  be  obtained  from  very  small  plots  In  poorly  reproduced 
reports,  and  hence  could  have  substantial  errors. 

In  the  following  chronological  review,  no  direct  quantitative  com¬ 
parison  of  any  two  of  these  experiments  will  be  attempted.  Such  an 
attempt  Is  precluded  by  the  variety  of  different  flow  parameters  and 
excitations  Involved.  Only  qualitative  consensus  on  the  response  to 
unsteadiness  will  be  sought. 

The  first  systematic  study  of  a  flat-plate  turbulent  boundary  layer 
subjected  to  uniform  sinusoidal  free-stream  oscillations  about  a  con¬ 
stant  mean  was  performed  by  Karlsson  [1959].  Measurements  of  tne  mean 
velocity,  the  amplitudes  of  the  In-  and  the  out-of-phase  components  of 
the  first  harmonic  of  the  periodic  velocity  and  the  intensity  of  higher 
harmonics  and  turbulence  were  made  using  a  hot-wire  anemometer  (HWA), 
over  a  wide  range  of  frequencies  and  amplitudes.  Karlsson  found  the 
mean  velocity  to  be  Independent  of  both  amplitude  and  frequency  of  exci¬ 
tation.  The  near-wall  velocity  always  led  the  free-stream  oscillations 
In  phase,  and  for  the  high-amplitude,  high-frequency  excitation,  there 
were  periodic  flow  reversals  near  the  wall. 

Hussain  and  Reynolds  [1970a,  1970b,  1972]  and  Reynolds  and  Hussain 
[1972]  reported  on  the  effect  of  a  weak  planar  disturbance  on  a  fully 
developed  channel  flow.  The  measurements  of  the  amplitude  and  the  rela¬ 
tive  phase  of  the  induced  streamwlse  velocity  components  were  made  in 
air  using  HWAs.  The  wave-induced  velocities  were  found  to  be  small  and 
to  decay  exponentially  In  the  mean-flow  direction.  The  decay  rates 
Increased  while  the  wave  lengths  decreased  as  the  frequency  was  In¬ 
creased.  The  authors  concluded  that,  even  for  such  small  disturbances, 
the  Induced  oscillations  In  the  Reynolds  stresses  cannot  be  neglected 
for  predicting  the  observed  wave  behavior. 

Acharya  and  Reynolds  [1975]  subjected  the  same  fully  developed 
channel  flow  to  periodic  variations  in  the  through- flow  rate.  They 
found  that  the  time-averaged  velocity  and  the  longitudinal  turbulence 
Intensity  were  practically  the  same  as  for  the  corresponding  steady 


flow.  The  amplitude  and  the  phase  of  the  periodic  component  of  the 
velocity  were  constant  over  the  central  portion  of  the  channel.  Indi¬ 
cating  a  "slug-like"  response.  All  the  variations  were  confined  to  a 
thin  region  near  the  wall,  and  the  behavior  for  the  two  frequencies 
studied  was  very  different.  This  was  attributed  to  the  fact  that  the 
higher  frequency  was  very  close  to  the  turbulent  "burst"  frequency.  It 
must  be  noted  that  the  very  small  signal-to-noise  ratios  resulting  from 
the  small  amplitudes  of  excitation  employed  caused  considerable  scatter 
in  the  near-wall  data.  Several  of  the  common  turbulence  closure  schemes 
were  found  to  perform  poorly  in  post-dieting  the  flow,  and  the  authors 
conjectured  that  higher  levels  of  closure  might  be  needed  to  Improve  the 
computations. 

Schachenmann  and  Rockwel 1  [1976]  reported  measurements  on  a  mean 
adverse  pressure  gradient  boundary  layer  In  a  conical  diffuser  subjected 
to  sinusoidal  variations  in  core  velocity.  Over  the  range  of  frequen¬ 
cies  studied,  they  found  no  detectable  effect  o?  unsteadiness  on  the 
mean  velocity  profiles.  They  found  that  the  amplitude  of  the  velocity 
oscillations  In  the  boundary  layer  could  exceed  the  core  velocity  ampli¬ 
tude  by  as  much  as  100%.  For  low-frequency  excitation,  the  phase  of 
velocity  In  the  boundary  layer  was  consistently  ahead  of  the  core  veloc¬ 
ity  phase;  for  high-frequency  excitation,  the  boundary- layer  velocity 
phase  lagged  that  of  the  core  flow,  and  the  amplitude  could  exhibit  more 
than  one  peak.  It  must  be  observed,  however,  that  there  are  only  a  few 
points  in  the  profiles,  and  the  data  scatter  is  quite  large  for  ampli¬ 
tudes  and  phases. 

The  group  at  ONERA/CERT  In  France  has  been  very  active  In  studies 
on  a  turbulent  boundary  layer  subjected  to  sinusoidal  free-stream  oscil¬ 
lations.  Coustelx,  Desopper,  and  Houdevllle  [1977],  Houdeville  and 
Coustelx  [1979],  and  Coustelx,  Houdevllle,  and  Javelle  [1981]  reported 
extensive  measurements  of  phase-averaged  velocity  and  turbulence  normal 
and  shear  stresses  using  HWAs  and  a  Laser  Doppler  Anemometer  (LDA)  on 
turbulent  boundary  layers  with  and  without  mean  adverse  pressure  gradi¬ 
ents.  The  mean- flow  development  was  unaffected  by  unsteadiness.  The 
phase  and  the  amplitude  of  boundary  layer  parameters,  like  the  displace¬ 
ment  and  the  momentum  thicknesses,  exhibited  an  oscillatory  behavior 


with  respect  to  the  Strouhal  number  based  on  the  length  of  boundary 
layer  development.  Hie  phase- averaged  boundary  layer  profiles  followed 
the  logarithmic  law  of  the  wall;  and  the  wall  shear  velocity,  computed 
using  Clauser  technique,  showed  Increase  In  amplitude  while  its  phase 
reduced  to  zero  as  the  Strouhal  number  increased.  The  velocity  phase 
plots  showed  a  local  maximum  near  the  wall  but  varied  considerably  very 
close  to  the  wall.  The  authors  conjectured  that  a  universal  near-wall 
law,  such  as  for  steady  flows,  was  perhaps  not  valid  for  unsteady  flows. 

Patel  [1977]  reported  Investigations  on  a  flat-plate  turbulent 
boundary  layer  subjected  to  harmonic  oscillations  of  the  free  stream 
velocity  associated  with  a  traveling  wave.  The  mean  flow  and  turbulent 
Intensity  were  unaffected  by  the  unsteadiness.  The  amplitude- ratio 
distributions  across  the  boundary  layer  were  found  to  be  independent  of 
the  free-stream  velocity  amplitude.  For  low-amplitude  cases,  this  was 
true  of  the  phase-angle  profiles  as  well.  For  higher- amplitude  cases, 
the  phase-lag  levels  were  found  to  Increase  with  the  free-stream 
velocity  amplitude.  A  linearized  high-frequency  analysis  using  an  eddy- 
viscosity  turbulence  model  compared  reasonably  well  with  the  experi¬ 
ments.  More  Important,  the  analysis  showed  the  dominant  Influence  of 
the  traveling-wave  convection  velocity  on  the  response  of  the  turbulent 
boundary  layer  to  free-stream  velocity  oscillations. 

Following  Patel's  work,  ken i son  [1977]  studied  the  effect  of  a 
traveling  wave  on  the  behavior  of  a  separating  turbulent  boundary  layer. 
The  mean  distributions  of  velocity  and  turbulence  profiles  and  skin 
friction  were  unaffected  by  the  unsteadiness  for  the  range  of  frequen¬ 
cies  and  amplitudes  studied.  The  velocity  amplitudes  In  the  boundary 
layer  were  dependent  on  the  frequency  and  the  distance  along  the  plate. 
There  was  a  thin  region  of  periodic  reversed  flow  ahead  of  the  mean 
separation  point,  but  It  had  no  significant  effect  on  the  variation  of 
momentum  thickness  or  shape  factor.  The  phase  of  the  velocity  oscilla¬ 
tions  was  always  a  lag  with  respect  to  that  of  the  free  stream.  The  lag 
levels  generally  Increased  with  downstream  distance  and  were  much  larger 
than  the  results  from  Patel  In  the  same  facility,  but  with  no  streamwlse 
pressure  gradient. 
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Ramaprlan  and  Tu  [1979,  1980,  1981]  have  reported  experiments  on  a 
fully-developed  pipe  flow  with  sinusoidal  variation  In  through-flow 
rate.  Measurements  were  made  In  oil  and  water  using  U)A,  and  a  range  of 
frequencies  was  covered.  The  authors  found  that  at  the  highest  fre¬ 
quency  of  oscillation  the  Reynolds  stresses  remained  frozen  at  an  aver¬ 
age  value  over  the  oscillation  cycle.  The  time-averaged  mean  velocity 
exhibited  a  point  of  Inflection  and  an  Increase  In  time-mean  shear 
stress  at  the  wall.  The  authors  argued  that  the  frequency  of  oscilla¬ 
tion  was  very  close  to  the  turbulent  burst  frequency,  causing  the  tur¬ 
bulence  to  Interact  with  the  organized  fluctuations,  resulting  In  a 
change  In  the  mean-velocity  profile.  It  must  be  noted,  however,  that 
this  result  was  based  on  a  single  data  point  near  the  wall,  and  the 
time-mean  shear  stress  at  the  wall  was  found  by  the  gradient  of  the  mean 
velocity  profile  at  the  wall,  a  highly .uncertain  procedure. 

Binder  and  Kuney  [1981]  studied  the  near-wall  region  of  a  fully- 
developed  channel  flow  subjected  to  small  periodic  velocity  oscilla¬ 
tions,  using  LDA.  The  time-averaged  velocity  and  turbulent  Intensity 
were  unaffected  by  the  oscillations.  For  high  frequencies,  the  near¬ 
wall  amplitude  and  the  phase  were  found  to  follow  the  laminar  Stokes 
solution. 

Shemer  [1981]  reported  measurements  on  a  pulsating,  fully-developed 
pipe  flow  In  both  laminar  and  turbulent  regimes  at  the  same  Reynolds 
number.  For  both  regimes,  he  found  that  the  time-mean  properties  of  the 
flow  were  Independent  of  the  oscillations.  The  radial  dependence  of  the 
amplitude  and  the  phase  of  the  velocity  oscillations  were  found  to  be 
different  for  laminar  and  turbulent  flows.  The  maximum  of  the  velocity 
amplitude  for  the  turbulent  flow  occurred  at  the  center  for  low  frequen¬ 
cies,  but  shifted  toward  the  wall  as  frequency  Increased.  The  near-wall 
velocity  oscillations  always  lagged  that  at  the  center  of  the  pipe,  and 
the  lag  Increased  with  frequency. 

As  an  extension  of  their  experiments  on  a  steady,  separating,  tur¬ 
bulent  boundary  layer,  Simpson  [1977]  and  later  on  Simpson  et  al.  [1980, 
1983]  studied  the  effect  of  a  low-frequency  sinusoidal  variation  of  the 
free-stream  velocity  on  the  same  flow.  Extensive  measurements  were  made 
of  mean  and  turbulent  quantities  and  the  forward-flow  fraction  using  a 
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number  of  measurements.  Including  HWA,  LDA,  and  thermal  tufts.  Upstream 
of  detachment,  the  mean  profiles  were  unaffected  by  the  unsteadiness, 
and  the  turbulence  structure  was  said  to  be  quasi-steady.  During  the 
detachment  process,  large  amplitude  and  phase  variations  developed 
throughout  the  flow.  The  anomalous  near-wall  phase  variations  In  the 
upstream  flow  were  attributed  to  small  oscillations  In  the  probe  posi¬ 
tion  with  respect  to  the  wall  at  the  exciting  frequency. 

Finally,  some  review  articles  of  Interest  should  be  mentioned.  For 
an  excellent  overview  of  the  many  different  areas  of  unsteady  fluid 
flows  and  applications,  see  McCroskey  [1977].  Tellonls  [1977]  and 
Coustelx  [1979]  have  reviewed  both  experiments  and  calculations  of 
unsteady  boundary  layers,  both  laminar  and  turbulent.  A  comprehensive 
review  of  existing  unsteady  turbulent  boundary  layer  experiments  is 
given  by  Carr  [1981a,  1981b].  The  compilation  Included  descriptions  of 
experimental  apparatus,  flow  conditions,  summaries  of  acquired  data,  and 
significant  conclusions.  Carr  has  created  a  magnetic  data-tape  library 
containing  the  data  from  these  experiments.  The  above  review  and  the 
data  library  should  prove  to  be  Invaluable  for  present  and  future  un¬ 
steady,  turbulent  boundary  layer  flow  experimenters  and  computors.  The 
present  data  are  being  added  to  Carr's  library. 

1.3  Summary  of  Previous  Research  and  Unanswered  Questions 

The  only  clear  conclusion  to  be  drawn  from  the  above  works  is  that, 
over  a  fairly  large  range  of  amplitudes  and  frequencies  of  excitation, 
the  time-averaged  flow  Is  practically  unaffected  by  the  Imposed  organ¬ 
ized  unsteadiness.  Several  unanswered  questions  remain. 

•  What  are  the  parameters  that  govern  unsteady  flows? 

Definitely,  a  frequency  parameter  In  the  form  of  a  Strouhal  number  and 
an  amplitude  parameter  should  be  on  the  list.  Depending  on  the  choice 
of  the  length  scale  (for  example,  streamwlse  distance,  boundary  layer 
thickness,  boundary  layer  Integral  thicknesses  are  all  possible  candi¬ 
dates),  different  Strouhal  numbers  can  be  defined.  None  seems  to  emerge 
as  a  clear  choice. 


•  What  Is  the  Interaction,  If  any,  between  the  organized 
fluctuations  and  the  disorganized  turbulent  fluctuations? 


No  clear  answer  to  this  very  important  question  emerges  from  the  avail¬ 
able  literature.  There  is  a  more  important  and  related  question:  how 
would  one  go  about  detecting  any  such  interaction?  It  would  seem  that 
any  correlations  between  an  organized  quantity  and  a  turbulence  quantity 
that  one  might  look  for  to  detect  interactions  would  be  suppressed  by 
the  phase-averaging  process  employed  by  almost  all  researchers. 

•  At  what  range  of  frequencies  can  one  expect  interactions 
between  the  organized  and  the  turbulent  fluctuations? 

There  is  an  expectation  among  most  researchers  that  the  frequencies  of 
excitation  have  to  be  in  the  range  characteristic  of  turbulent  fluctua¬ 
tions  before  any  interactions  will  take  place.  The  Strouhal  number 
based  on  the  boundary  layer  thickness  for  the  turbulent  "burst"  fre¬ 
quency  is  about  1.26,  based  on  the  results  from  Rao,  Narasimha,  and 
Narayanan  [1971].  With  a  few  exceptions,  in  the  experiments  above,  the 
excitation  frequencies  were  considerably  lower  than  the  "burst"  fre¬ 
quency.  No  definite  reasons  or  arguments  emerge  from  the  literature, 
explaining  why  sub-harmonic  interactions  are  impossible  or  unlikely. 

The  works  reviewed  above  point  to  the  lack  of  sufficient  experi¬ 
mental  data  on  unsteady  turbulent  boundary  layers.  Some  important 
considerations  in  meeting  these  data  requirements  are  discussed  next. 


1.4  Data  Requirements  in  Unsteady,  Turbulent  Boundary  Layers 

Two  needs  emerge  from  a  careful  review  of  the  above  literature: 

•  need  for  very  carefully  done  experiments  on  very  clearly 
defined  unsteady  turbulent  flows; 

•  need  for  systematic  and  careful  evaluation  of  the  perfor¬ 
mance  of  existing  turbulence  models  in  such  flows. 

Some  Important  considerations  in  fulfilling  the  first  of  these  needs  in 
turbulent  boundary  layers  is  addressed  below. 

A  good  physical  understanding  of  unsteady,  turbulent  boundary  lay¬ 
ers  can  be  gained  only  through  careful  study  of  their  dynamic  behavior 
when  subjected  to  different  types  of  unsteadiness.  The  organized 
unsteadiness  Imposed  on  a  turbulent  flow  can  be  of  several  different 
types.  It  could  be  periodic  or  transient.  If  periodic,  then  it  could 


be  of  different  wave  forms.  The  wave  forms  could  correspond  to  a 
standing  wave  pattern  or  could  be  traveling  at  a  speed  more  or  less  than 
the  reference  flow  speed.  If  the  phenomenon  In  question  were  a  linear 
one,  one  could  construct  the  response  for  a  given  type  of  unsteadiness 
from  a  superposition  of  the  responses  to  elementary  harmonic  excitations 
over  the  relevant  frequency  range.  If  the  phenomenon  Is  a  nonlinear 
one,  the  response  could  be  potentially  very  different  for  different 
types  of  excitations.  Turbulence  Is  essentially  a  nonlinear  phenome- 
non.  There  Is  no  reason  to  expect  that  the  response  of  turbulent  flows 
to  organized  unsteadiness  will  be  linear  except  perhaps  for  very  small 
amplitudes.  Viewed  from  this  perspective,  the  available  data  on  un¬ 
steady  turbulent  boundary  layers  are  very  meager. 

For  a  given  type  of  Imposed  unsteadiness,  it  is  important  to  cover 
a  large  region  in  the  parameter  space  to  identify  different  flow  re¬ 
gimes,  if  any.  In  most  of  the  experiments  reviewed  above,  independent 
variation  of  parameters  was  difficult.  If  not  impossible.  In  particu¬ 
lar,  most  of  the  above  experiments  do  not  cover  a  wide  range  in  the  fre¬ 
quency  parameter.  Even  for  those  that  do  (Karlsson  [1959]  and  Schachen- 
mann  and  Rockwell  [1976]),  the  data  are  available  only  in  the  form  of  a 
few  plots  In  the  cited  references;  moreover,  these  results  suffer  from 
lack  of  definition  and  deficiencies  in  instrumentation  and  data- 
processlng. 

In  developing  an  understanding  of  the  dynamics  of  unsteady  turbu¬ 
lent  boundary  layers.  It  Is  very  Important  to  minimize  effects  due  to 
other  causes.  All  of  the  works  reviewed  above  are  characterized  by 
unsteadiness  throughout  the  flow  field.  The  effects  of  organized  un¬ 
steadiness  on  transition  to  turbulence  and  the  subsequent  development  of 
the  turbulent  flow  are  wholly  untested,  and  remain  open  questions.  Thus 
It  Is  not  clear  how  much  of  the  effects  seen  In  the  above  works  should 
be  attributed  to  unsteady  transition.  Certainly,  from  the  computational 
point  of  view.  It  would  be  desirable  to  start  from  steady  conditions, 
and  then  Impose  unsteadiness. 

The  present  experiments  were  specifically  designed  to  address  many 
of  the  above  requirements. 
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1.5  Objectives 

The  primary  motivations  for  the  Stanford  Unsteady  Turbulent  Bound¬ 
ary  Layer  Research  Program  are  two- fold: 

•  to  explore  the  possibilities  of  controlling,  reducing,  or 
eliminating  unsteady  fluid  dynamic  effects  such  as  dy¬ 
namic  stall  by  carefully  tailored  and  externally  imposed 
excitations,  and 

•  to  provide  design  guidelines  for  the  improvement  in  per¬ 
formance  of  fluid  dynamic  devices  with  inherent  unstead¬ 
iness,  such  as  turbomachinery  and  internal  combustion 
engines. 

The  overall  objectives  of  the  Stanford  Program  are  to  provide  a 
fundamental  understanding  of  the  fluid  dynamical  and  heat-transfer 
characteristics  of  unsteady,  turbulent  boundary  layers  in  general,  and 
to  acquire  a  comprehensive  data  base.  From  the  overall  objectives  of 
the  Program,  the  following  specific  objectives  of  the  research  reported 
here  were  formulated: 

•  To  understand  the  characteristics  of  a  turbulent  boundary 
layer  with  well-defined  initial  and  boundary  conditions. 

•  To  provide  a  definitive  data  set  on  such  an  unsteady, 
turbulent  boundary  layer  that  will  serve  as  a  test  case 
for  evaluating  the  performance  of  calculation  procedures 
and  as  Input  for  turbulence  models  of  such  flows. 

A  brief  outline  of  how  these  objectives  were  met  is  given  next. 

1.6  Sumnary  of  Present  Work 

To  meet  the  above  objectives,  a  closed-loop  water  channel  was 
specially  constructed  for  the  experiments.  In  the  upstream  development 
section  of  the  channel,  a  turbulent  boundary  layer  was  grown  under 
steady,  zero-pressure  gradient  conditions.  This  provided  a  well-defined 
Inlet  condition  to  the  test  section,  where  the  boundary  layer  was  sub¬ 
jected  to  a  sinusoidally  varying,  adverse,  free-stream  gradient.  A 
laboratory  computer  was  set  up  to  automate  the  data  acquisition  and 
control  procedure. 

Two  major  series  of  experiments  with  relative  amplitudes  (of  free- 
stream  oscillations  at  the  end  of  the  test  section)  of  5%  and  25%  were 
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carried  out.  Several  frequencies  and  X  stations  were  covered  In  each, 
to  provide  a  wide  range  of  Strouhal  number.  Measurements  of  phase- 
averaged  streamwlse  velocity  and  square  of  streamwlse  turbulent 
fluctuations  were  made,  from  which  other  quantities  of  Interest  were 
computed. 

The  time-averaged  flow  was  found  to  be  Insensitive  to  the  Imposed 
unsteadiness.  The  periodic  flow  was  strongly  affected  by  frequency  and 
showed  two  distinct  regimes  of  behavior.  The  high-frequency  results 
agreed  well  with  the  Stokes  solution.  .Over  the  range  of  frequencies  and 
amplitudes  studied,  no  dramatic  Interactions  between  the  organized 
oscillations  and  the  disorganized  turbulent  fluctuations  were  found. 

An  extensive  data  set  has  been  obtained.  This  Is  documented  In 
this  report  and  a  separate  microfiche  report  LJayaraman,  Parikh  and 
Reynolds,  1982 j.  A  digital  magnetic  tape  containing  all  the  data  Is 
available  from  Dr.  Lawrence  W.  Carr,  U.  S.  Ann/  Aeromechanics  Labora¬ 
tory,  Anes  Research  Center,  MS  215-1,  Moffett  Field,  California  94U35. 


Chapter  2 

THEORETICAL  BACKGROUND 


2.1  Triple  Decomposition 

Following  Hussain  and  Reynolds  [1970],  any  flow  variable  f(x_,t) 
such  as  velocity  or  pressure,  can  be  decomposed  into  three  components: 

f(x,t)  =  f(x)  +  f(x,t)  +  f'(x,t)  (2.1) 

Here  f(_x)  is  the  time-averaged  or  mean  value  of  f(x^t)  at  a  given 
location,  f(x,t)  is  the  organized  response  component  due  to  the  im¬ 
posed  organized  unsteadiness,  and  f'(x_,t)  is  the  turbulent  fluctua¬ 
tion.  Several  unsteady  turbulent  flow  Investigations  to  date  have  used 
the  above  decomposition. 

2.2  Averaging 

In  order  to  separate  the  three  components  of  a  flow  variable,  two 
different  averaging  procedures  are  employed,  namely,  time  average  and 
ensemble  average.  The  time  average  is  conventionally  defined  as  fol¬ 
lows: 

T 

f (_x)  *  lim-i  f  f (_x ,t)  dt  (2*2) 

T+-  T  J 
o 

When  the  variable  is  sampled  at  discrete  times,  the  time  average  is  the 
average  of  a  large  number  of  samples  in  a  record  of  sufficiently  long 
duration.  In  other  words, 

N 

f(x)  «  llm^^ftx,^)  (2.3) 

^  1*1 

where  t^  are  the  discrete  times  (usually  equally  spaced)  when  samples 
were  taken.  It  Is  also  understood  that  the  total  time  duration  In  which 
all  the  samples  were  taken  Is  large  compared  to  both  the  Imposed  excita¬ 
tion  period  and  the  longest  turbulent  fluctuation  period. 

In  order  to  extract  the  organized  component,  an  ensemble-averaging 
procedure  is  used.  An  ensemble  average  is  defined  as  the  average  over  a 


large  number  of  realizations  of  the  same  phenomenon.  This  Is  expressed 
as: 

N 

<  f(x,t)  >  -  llm-jj^jf^x.t)  (2.4) 

^  1*1 

where  1  Identifies  a  particular  realization.  Further,  each  realiza¬ 
tion  Is  understood  to  have  begun  at  t  *  0. 

For  an  Imposed  periodic  excitation,  the  organized  response  compo¬ 
nent  will  also  be  periodic  with  the  same  period;  thus  the  realization, 
1,  will  occur  at  the  same  point  In  each  cycle,  and  the  organized 
component  can  be  called  the  periodic  component.  The  ensemble  average 
for  this  case  becomes 

N-l 

<  f(x,t)  >  *  11m  4  V"  f(x,t  +  n-c)  (2.5) 

N~  "  “  " 

n*0 

where  x  Is  the  time  period  of  the  Imposed  excitation. 

It  Is  easier  to  think  of  the  ensemble  averaging,  for  periodic  exci¬ 
tations,  as  phase  averaging.  A  phase  average  Is  defined  as  the  average 
of  samples  from  a  large  number  of  cycles  at  a  fixed  cycle  phase  angle. 
That  Is, 

N-l 

<  f(x,0)  >  -  llm-J^  f(x,0  +  2*  n)  ,  0<e<2*  12.6) 

N-mo 

n*0 

where  e  Is  the  cycle  phase  angle.  Definitions  (2.5)  and  (2.6)  are 
entirely  equivalent.  In  practice,  N  is  always  finite  but  large. 

By  definition. 


<  f(x,t)  >  *  f (x)  +  f(x,t)  (2.7) 

and  the  time  average  (or,  equivalently,  the  average  over  all  phases)  of 
the  phase  average  Is  the  mean.  That  is, 

<'  'f(x.'t)'  >■  -  T(x)  (2.8) 


Therefore 


f(x,t)  -  <  f(x,t)  >  -  f(x)  (2.9) 

and 

f'(x.t)  »  f(x,t)  -  <  f(x,t)  >  (2.10) 

From  the  above  definitions,  the  following  useful  relations  can  be  devel¬ 
oped: 

<  f  g  >  =  f  <  g  > 

<?>  «  7 

<7g>  =  7<g>  (2.11) 

fg*  -  0 


<  f  >  *  0 
f  g  =  f  g 

f  *  0 

f'  -  0 


<  fg'  >  -  0 


2.3  Fourier  Analysis 

Fourier  analysis  Is  useful  In  determining  the  amplitudes  and  the 
phase  angles  of  the  harmonic  components  of  the  periodic  response.  Any 
periodic  waveform  can  be  split  up  into  its  Fourier  components  as  fol¬ 
lows: 

m 

?(x,t)  ■  ^2  An,fM  c°s  |^nwt  +  ^2>12^ 

n»l 

where  A^fJxJ  and  +n  f(x)  are  the  amplitude  and  the  phase  of  the 
nth  harmonic  of  the  periodic  component  f(>c,t)  at  that  location  and 

u  -  ~  (2.13) 

T 

where  t  is  the  time  period  of  the  Imposed  periodic  excitation. 
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Complete  description  of  f  requires  specification  of  Ap  f  and 
for  all  values  of  n.  Only  the  first  few  harmonics  are  usually 
Important  In  the  response  of  physical  systems  to  a  purely  sinusoidal 
excitation. 


2.4  Governing  Equations 

The  equations  of  motion  provide  considerable  Insight  Into  the 
physics  of  the  flow.  The  governing  equations  (using  Cartesian  tensor 
notation)  for  an  Incompressible  flow  of  a  Newtonian  fluid  are  as 
follows: 


Continuity: 

Momentum: 

Let 


du^ 

ax. 


du^  bU| 

-  +  u  .  - 

at  j  dx . 

j 


a2u. 

i  &p  ^  _ i_ 

p  ax.  v  bx.&x. 
1  j  j 


(*.14) 


12.1b) 


u  ■  u  +  u  +  u1 
1  1  1  1 


p  ♦  p  +  p 


(2.16a) 

(2.16b) 


Substituting  Eqs.  (2.16)  Into  (2.14)  and  performing  phase  and  time  aver¬ 
aging. 


8ui 

6xi 

0 

(2.17a) 

6“i 

8xi 

0 

(2.17b) 

au; 

axi 

0 

(2.17c) 

Substituting  Eqs.  (2.16)  Into  (2. IS), 


it  (ui +  ui +  u;> +  irL(uj +  uj  +  uj)(ui +  ui +  vj 

J 

2 

"  *  p*  p,)  *  ’inr  (ui  *  “i  *  “i> 

1  J  J 

or 

—  (i  +  u '  +  ui)  +  — —  (UM.  +  u.u.  +  u.u'.  +  u.u.  +  u  u.  +  u  u'. 

at  '  1  1  r  5x.  1  1  j  1  j  1  j  1  j  1  j  1  j 

J 

*  ui"j +  ui“i  *  “i"?  ■  'Isq-  lp * p,)  * 


*  'HIT  <U1  *  “i  +  “S’ 

J  J 


Phase  averaging  (2.18), 


— 1  +  (u.u  ♦  u.u  +  u  u  +  u  u*  +  <  uV.  >) 

at  ax  1  j  1  j  i  j  i  j  i  j 

J 

2 

la,—  a  ~ . 

■  -p«r(|,*|,1  +  v»rM:(uitui) 

1  J  J 


Time  averaging  (2.18), 

—  A, 

-JL-  (uu  +  uu  ♦u'u1)  “  -  —  ^£_  +  v - 

»x  lVj  Vj  uiV  p  &xi  axjaxj 

Subtracting  (2.20)  from  (2.19), 

«w  2~ 

5u  ,  -  a  u. 

at  axj  1  i  j  Vj  i  j  IT  pax1  ^axj 

Equation  (2.20)  can  be  expressed  in  the  more  usual  form  as 


a  . —  ,  l  ep  . 

■  ‘  pSX^ 

Equation  (2.21)  can  be  rewritten  as 


«  u.  _  _ 

- - —  (u'u*  +  u.u  .)  (. 

ax. ax.  ax.  v  i  j  i  j 

j  j  j 


8u.  „  .  ~  5  u  , -  , — ^ 

— 1  +  J_  (u  S'  +  u  u  )  -  -!*£_♦  v - -  -  _L  (u„u4  +  u'u1.) 

at  axj  '  i  j  i  j;  p  bx1  ax^Xj  ax^  v  i  j  i  jj 


For  a  2-D  phase-averaged  flow,  let 


1 

Equations  (2.17)  can  be  rewritten  as 


(u,v,w) 

(U.V.O) 

(x,y,z) 


bU  bV 
—  +  —  ■  o 

bx  by 

(2.24a) 

$£**1  -  o 

9x  9y 

(2.24b) 

Equation  (2.22)  can  be  expressed  as 


f;  (D  T)  ♦  (57) 


(U’U'J  -  -gy  (uV) 


and 


h lv  “>  ♦  h <y  v) 


Up,  A2U  „  .2u\  9 

Vj'* 

8  ,  ~  8  ~ 

*  9*  U  U  "  9y  U  V 

.!«£♦  ,/»!*. »_  (Tv)  .1- 

U2  9,7  “  •» 


(2.25a) 


8  ~  ~  8 
-  —  (v  u)  -  —  (v  v) 
8x  '  '  by 


(vV) 

(2.25b) 


Equation  (2.23)  can  be  rewritten  as 


1  bp  /  8Zu  8^ 

"  p  bx  +  Vl7T  +  — 2 
\bx  by 


5  ^  5  m*  **  1  1,1  ^ 

-  ~  (u  u  +  u'u' )  -  —  (u  v  +  uV) 
bx  '  '  by  v  ' 


(2.26a) 


and 


M  ^  2^> 
tv  8  ^  ^  ft  —  ^  1  ftp  i  ft  v  ft  v 

« * ?; (v “ +  * u) * 2 5y (< *>  *  *p»r  nr?  *77 

\»*  ay  / 

-  •j-  (u  v  ♦  uV)  -  -5-  (v  v  ♦  vV)  (2.26b) 
bx  '  by 
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2.5  Simplification  of  Governing  Equations 


The  purpose  of  this  section  is  to  simplify  the  governing  equations 
given  in  Section  2.4  for  the  special  case  of  a  2-D  Incompressible  bound¬ 
ary  layer.  This  Mill  be  achieved  by  performing  an  order  of  magnitude 
analysis  on  the  equations  and  eliminating  terms  that  are  small  compared 
to  others.  For  convenience,  the  simplified  equations  are  summarized  in 
Table  2.1.  The  details  of  the  analysis  are  given  In  Sections  2.5.1, 
2.5.2,  2.5.3  and  2.5.4. 

2.5.1  Primary  Estimates 

In  performing  the  order-of-magnltude  analysis,  variables  will  be 
replaced  by  "typical"  (order  of  magnitude)  values.  The  changes  in  vari¬ 
ables  are  expected  to  be  of  the  same  order  as  typical  values.  In  the 
following,  the  symbol  Is  used  to  Indicate  a  typical  value  of 

order-of-magnltude  accuracy  without  regard  to  the  sign,  this  should  not 
be  confused  with  the  "~"  used  over  a  periodic  variable.  The  symbol 
"A"  In  front  of  a  variable  Is  used  to  Indicate  changes  in  that  vari¬ 
able. 

For  a  2-D  Incompressible  boundary  layer  flow,  the  following  orders 
of  magnitude  are  assumed: 

Ax  ~  L  for  both  mean  and  periodic  variables. 

At  ~  1/u), 

Ay  ~  6  for  mean  variables, 

~  6*  for  periodic  quantities, 

U,  /IT  ~  ue, 

U,  AU  ~  A, 

where 

L  Is  a  typical  length  scale  In  the  streamwlse  direction  over 
which  typical  changes  In  mean  and  periodic  variables  take 
place. 


i 

i 


u  Is  the  circular  frequency  of  oscillation, 

6  Is  the  mean  boundary  layer  thickness  (6  «  L), 

6‘  Is  a  typical  length  scale  In  the  normal -to-wa  11  direction  over 
which  typical  changes  In  periodic  variables  occur  (normally  6* 
1  6), 

Ue  Is  the  mean  free- stream  velocity,  and 
A  Is  the  amplitude  of  the  free-stream  velocity  oscillations. 

It  Is  further  assumed  that 

A(u'u')  ,  A(uV)  ~  e2U2 

e 

and 

A(u'u')  ,  A(u'v‘)  ~  eZo  U2  , 

G 

where 

t2  is  a  multiplicative  factor  (usually  «  1)  whose  order  of 
magnitude  will  be  determined  later  In  the  analysis  and 

«  -  A/Ue. 

The  order  of  magnitude  of  the  changes  In  (u'u1)  and  (u'v1)  may  be 
arrived  at  by  the  following  reasoning.  Assume  that,  under  quasi-steady 
(u  ■  0)  conditions. 


2  2 
<  u'u1  >  ~  c  <  u  > 

e 

where  <  ue  >  Is  phase-averaged  free-stream  velocity  and 

<  u  >  -  U  +  A  cos(<rft) 
e  e 


Therefore , 


<v 


U2  *  ZAU  cos  (ut)  +  A2  COS2 
e  e 

U2  |~1  +  2«  cos  (wt)  +  a2  cos2 


(ut) 

(*>] 


18 


Usually  o  <  1  and  a2  Is  small  compared  to  a.  Thus 


(u'u')  ,  A(u'u')  ~  z  a  . 

A  similar  reasoning  holds  for  A(u'v').  Experimental  data  presented  in 
Chapter  4  show  that  the  estimate  for  A(uV)  is  good  over  all  the  fre¬ 
quencies  studied. 

Each  of  Equations  (2.24),  (2.25),  and  (2.26)  will  be  analyzed  indi¬ 
vidually  in  the  following  sections. 


2.5.2  Analysis  of  Continuity  Equations 
Consider  Eq.  (2.24a): 

au  Ue  .  a7  aV 

—  —  and  —  m  — 

ax  L  ay  6 

Since  both  terms  must  be  significant. 


AV  ~  1 


Consider  Eq.  (2.24b): 


au 

ax 


A 

I 


and 


This  implies  that 


Term  4  Is  negligible  compared  to  Term  5,  and  Term  6  Is  negligible  com¬ 
pared  to  Term  7,  since  ( 6/L)  «  1.  For  the  viscous  term  to  be  signifi¬ 
cant  (Term  5),  we  require  that 
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For  the  pressure  gradient  term  to  be  significant.  Me  note  that 


Further  note  that  Terms  8  and  9  are  of  order  o£  compared  to  Terms  1 
and  2  and  hence  are  negligible  even  for  a  as  great  as  0.3. 

Finally,  for  the  only  turbulent  term  to  be  of  any  significance 
(Term  7),  we  note  that 


2 

c 


Consider  Eq.  (2.25b);  there  are  nine  terms, 
magnitude  of  each  is  estimated  as  folloMS. 


Term  1: 


Term  2: 


U  fiU 
e  e 

■  ■  s 

LL 


,  6U  fill 
lee 


The  order  of 


Term  5: 


Term  8: 


Term  9: 


M<“  V) 


1  .  AS' 

ta  — 


A26* 


8  /~  *•» 
—  (v  v) 
8y 


1  6‘A  6'A 
6*  L  L 


.  %£» : 
2  2  6 
L  U 

U*6  A2  6* 

2  ~2  ~6~ 

L  U 

e 


With  the  above  estimates,  for  Eq.  (2.25b),  Terms  1,  2,  4,  5,  6,  8  and  9 
are  negligible  (for  a<  0.3)  compared  to  Term  7.  Thus 


8p 

ay 


pu 


or 


Therefore, 


<v?0>  -  «  1  • 


*£  .  0 
ay 


(2.27) 


and  Eq.  (2.25a)  simplifies  to 

-  -lJ.vO.|_(^)  ,.<0.3.  (2.28) 

This  equation  Is  the  same  as  the  stea<(y  boundary  layer  momentum  equa¬ 
tion. 


2.5.4  Analysis  of  Periodic  Momentum  Equations 

Equation  (2.26a)  has  11  terms.  Their  orders  of  magnitude  are 
estimated  as  follows: 


Term  1: 


Au 


Term  2: 
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S  «  1 


For  the  present  experiments,  L  ~  0.6  m,  Ue  ~  0.7  m;  the  above  re¬ 
quirement  translates  Into 

<d  «  1  ra d/s 

For  this  case,  the  unsteady  term  can  be  neglected  in  comparison  to  other 
terms.  Further,  if  viscosity  has  to  be  Important,  then  6‘  ~  6,  i.e., 
the  thickness  of  region  over  which  unsteady  effects  are  Important  is  of 

the  same  order  as  the  boundary  layer  thickness.  The  pressure  oscilla- 

2 

tlons  are  of  the  order  of  pAUe  *  pUga.  Thus  the  governing  equation  for 
u  becomes  (anticipating  the  result  for  p  to  be  shown  later): 


M  ">  * 


5  I 

w  (“  V)J 


(2.29) 


The  terms  in  square  brackets  are  of  order  a  and  hence  are 
probably  Important  only  for  a  >  0.1. 

Case  2  (very  high  frequencies) 


Ab>  » 


(i)L 

r  »  1 

e 

For  the  present  experiments,  this  means 


w  »  1  rad/s 

In  this  case,  for  the  viscous  term  to  be  significant, 

1  , 

*-  L  »'2  ~ 
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6'  1 
6-  ~  sl/2 

Clearly,  at  very  high  Strouhal  numbers,  the  unsteady  effects  are  con¬ 
fined  to  a  very  thin  region  near  the  wall. 

For  this  case  the  pressure  oscillations  are  of  the  order  of  pAuL  * 
pUgOS.  Thus  the  pressure  oscillations  increase  with  Strouhal  number. 
The  governing  equation  reduces  to  the  Stokes  equation,  as  shown  below: 


au 

at 


(2.30) 


Case  3  (Intermediate  frequencies): 


Aw  ~ 


wL  c  . 

\T  m  S  ~  1 
e 

For  this  case,  the  unsteady  term  has  to  be  added  to  Eq.  (2.29). 
~  5p  c**  . 

au  „  a  .  ~  ,  a  ~  ~  l  e  au  a  ,  ,  ,  v 

stt2«(uU)+»7(u,tuV)  '  -p»T+  V-T-*F(U,) 


-  L65(u  u)  +  iS  (u  V,J 


(2.31) 


Once  again,  the  terms  in  square  brackets  are  of  order  a  and  hence  need 
to  be  retained  only  for  a  >  0.1. 

„  To  complete  the  order  of  magnitude  analysis,  it  will  be  shown  that 
—  •  0  for  an  unsteady  boundary  layer.  Equation  (2.26b)  has  11  terms. 
Their  orders  of  magnitude  are  estimated  as  follows: 
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S  «  1,  6'  ~  6.  In  this  case. 


~  pAU  pU  a 
dp  e  e 


ay 


or 


Ap 


~  pH  a  r  «  1 

e  L 


Case  2: 


S  »  1;  (6*/6)  ~  1/S1/2  In  this  case. 


~  m 

dp  e  6 

ay  ~  P~T  61 


pAU 


ft 


or 


Ap 


•  *  t  ft) 


1/2 


pUg« 


since  usually 


U  L 
e 


»  S 


Case  3: 

S  ~  1;  6'  ~  6.  This  Is  the  same  as  Case  1 

Thus,  for  a  boundary  layer. 


Hence 


ap 

ay 


o 


ap 

ax 


as  had  been  used  earlier 


Chapter  3 

EXPERIMENTAL  SETUP  AND  PROCEDURES 


The  experimental  setup  consists  of  the  apparatus  and  the  data- 
acqulsltlon  and  control  system.  The  experimental  procedures  consist  of 
the  data-acqulsltlon  and  control  procedures  and  the  data- reduction 
procedures.  Each  one  of  these  will  be  briefly  described  below. 

3.1  Apparatus 

The  apparatus  consists  primarily  of  the  water  channel  and  the  os¬ 
cillation  mechanism  used  to  produce  the  free-stream  oscillations.  The 
objectives  behind  the  design  of  the  apparatus  will  be  described  first. 
After  a  description  of  the  water  channel  and  the  means  of  producing 
controlled  oscillations,  the  treatment  and  control  of  water  quality  and 
the  LDA  mount  will  be  briefly  touched  upon. 

3.1.1  Design  Objective 


Figure  3.1  shows  the  desired  free-stream  variation  In  the  experi¬ 
mental  apparatus.  In  the  upstream  development  section,  the  free-stream 
velocity  Is  steady  In  time  and  uniform  In  the  X  direction.  In  the 
test  section,  the  free-stream  velocity  decreases  linearly  with  X,  with 
the  rate  of  decrease  dependent  on  the  cycle-phase  angle.  At  0°  cycle- 
phase  angle,  the  free-stream  gradient  Is  0  and  Its  value  Is  the  same 
as  upstream.  In  the  recovery  section,  the  free-stream  velocity  is  again 
uniform  In  X  but  Is  a  function  of  cycle-phase  angle. 

During  operation,  the  free-stream  velocity  gradient  In  the  test 
section  oscillates  sinusoidally  at  the  desired  frequency,  between  0 
and  the  desired  maximum  value.  The  time-averaged  free-stream  velocity 
decreases  linearly  with  X.  The  amplitude  of  the  free-stream  oscilla¬ 
tions  Increases  linearly  with  X.  Thus  the  boundary  layer  In  the  up¬ 
stream  development  section  grows  under  zero  pressure-gradient  conditions 
as  a  steady,  flat-plate  boundary  layer.  In  the  test  section.  It  Is 
subjected  to  an  oscillatory  adverse  free-stream  gradient  with  a  nonzero 
mean. 
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The  experimental  parameters  that  can  be  varied  are  the  frequency  of 
oscillation  and  the  amplitude  of  free-stream  oscillations  at  the  end  of 
the  test  section.  The  actual  details  of  the  Implementation  of  the  de¬ 
sign  objective  In  the  test  apparatus  are  given  below. 

3.1.2  Stanford  Unsteady  Turbulent-8oundary-Layer-Research  Water 
Channel 

A  special  and  unique  water  tunnel  facility  in  which  unsteady  tur¬ 
bulent  boundary  layers  can  be  studied  In  great  detail  was  constructed 
for  the  research  program.  Water  was  chosen  as  the  working  fluid  for  the 
following  Important  reasons: 

•  The  fluid  velocities  required  are  low,  and  hence  the 
turbulence  frequencies  of  Interest  will  be  low.  Free- 
stream  oscillations  at  such  frequencies  (0-2  Hz)  can  be 
achieved  easily  with  relatively  simple  mechanical 
arrangements. 

•  Forward-scatter  Laser-Doppler  Anemometry  with  frequency¬ 
tracking  Is  easy  to  use  In  water. 

•  Water  permits  easy  visualization  of  the  flow  field. 

•  Unsteady  compressibility  effects  are  much  less  severe  In 
water  than  In  air. 

Figure  3.2  shows  a  schematic  diagram  of  the  Stanford  Unsteady 
Turbulent-Boundary- Layer-Research  water  channel.  Water  pumped  from  the 
sump  Is  supplied  to  the  overhead  tank  (section  1.5  x  l  m)  through  a 
standing  pipe  with  several  holes.  The  water  level  Is  maintained  con¬ 
stant  to  within  2  mm  ^0  by  an  overflow  weir.  An  open  rectangular 
header  (not  shown  In  the  figure)  Is  placed  Inside  the  tank  around  the 
entrance  to  the  downcomer.  The  header  has  a  honeycomb  for  one  of  Its 
vertical  faces  through  which  water  can  enter  the  header.  The  purpose  is 
to  make  the  flow  Inside  the  header  as  quiescent  as  possible. 

The  entrance  to  the  downcomer  Is  through  a  honeycomb.  At  the 

bottom  of  the  downcomer,  water  Is  turned  through  90°  by  a  triangular 
header  and  honeycomb  combination  Into  the  screen  box.  The  screen  box  Is 
made  of  plexiglas,  and  the  screens  are  made  of  stainless  steel  wire  (24  « 

gauge,  70%  porosity)  and  are  stretched  onto  plexiglas  frames.  There  are 
three  screens,  spaced  12  cm  apart,  and  they  are  Inserted  Into  the  box 
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much  like  printed  circuit  cards.  The  sides  of  the  box  can  be  disassem¬ 
bled  for  removal  and  cleaning  of  the  screens.  Unce  in  place,  the 
screens  can  be  put  into  tension  by  bolts  that  reach  from  the  sides  of 
the  box.  It  was  found  necessary  to  keep  the  screens  clean  and  their 
surface  as  planar  as  possible  to  achieve  very  good  spanwise  uniformity 
in  the  development  and  test  sections  (see  Figs.  4.1).  The  cross-section 
of  the  screen  box  is  1  m  x  1  m.  The  screen  box  is  followed  by  a  nozzle 
that  contracts  the  flow  in  both  directions  to  the  channel  cross- 
section  (35  cm  x  15  cm). 

The  nozzle,  the  overhead  tank,  the  downcomer,  and  tne  yo°  turn 
header  were  all  constructed  of  32  mm  steel  plates  reinforced  by  structu¬ 
ral  channels.  The  inside  surfaces  of  these  components  were  sandblasted, 
washed  with  acetone,  and  coated  with  epoxy  resin  to  prevent  rusting. 

The  nozzle  is  followed  by  a  2  m  long  development  section  of  cross- 
section  3b  cm  x  15  cm.  At  tne  end  of  the  nozzle  (i.e.,  at  the  begin¬ 
ning  of  the  development  section),  two  rows  of  closely  spaced  bleed  holes 
are  provided  for  removal  of  side-  and  top-wall  boundary  layers  from  the 
nozzle.  This  permits  growth  of  a  new  test  boundary  layer  on  the  top 
wall  of  the  development  section.  As  the  boundary  layers  grow,  the 
free-stream  velocity  in  the  core  increases  due  to  the  increasing  block¬ 
age.  Part  of  the  flow  needs  to  be  bled  off  to  keep  the  free-stream  vel¬ 
ocity  uniform.  For  this  purpose,  rows  of  bleed  holes  are  provided  in 
the  bottom  wall  every  15  cm  throughout  the  development  section.  The 
flow  through  each  bleed  hole  row  is  controlled  by  a  valve.  These  valves 
are  carefully  adjusted  to  achieve  a  very  uniform  free-stream  velocity  in 
the  development  section  (see  Fig.  4.2). 

The  top-  and  the  side-wall  boundary  layers  are  tripped  by  a  sand¬ 
paper-type  rough  tape  placed  12  cm  downstream  of  the  nozzle  exit.  At 
the  end  of  the  development  section,  the  side-wall  boundary  layers  are 
removed  by  splitter  plates,  and  the  channel  width  is  reduced  from  35  cm 
to  30  cm.  This  forces  growth  of  new  side-wall  boundary  layers  in  the 
test  section  that  will  be  thinner  than  the  top-wall  study  boundary  layer 
and  hence  be  more  resistant  to  separation.  The  amount  of  flow  removed 
Is  carefully  adjusted  to  give  a  smooth  flow  at  the  leading  edge  of  the 
splitter  plates,  visualized  by  dye,  injected  just  upstream  of  the 
piates. 


The  test  section  that  follows  the  development  section  is  61  cm  long 
and  30  cm  x  15  cm  in  cross  section.  The  bottom  wall  of  the  test  sec¬ 
tion  consists  of  a  porous  plate,  through  which  flow  can  be  extracted. 
The  test  section  Is  followed  by  a  short  recovery  section  and  then  an 
exit  section.  The  exit  section  also  has  a  porous  plate  similar  to  that 
In  the  test  section.  These  porous  plates  have  several  holes  spaced  1.3 
cm  apart  in  rows  spaced  every  2  cm  along  the  plates. 

To  have  uniform  suction  in  the  bottom-wall  porous  plate  in  the  test 
section,  the  pressure  drop  across  the  plate  should  be  large  compared  to 
the  change  in  static  pressure  in  the  channel  along  the  plate.  The  size 
of  the  holes  in  the  plates  is  chosen  for  each  amplitude  such  that  the 
above  requirement  is  met. 

The  top  wall  of  the  channel  forms  a  smooth  working  surface  for  the 
test  boundary  layer.  Probe  access  ports,  approximately  2  cm  in  diam¬ 
eter,  are  provided  along  the  centerline  of  the  channel  and  spanwise  at 
several  locations.  These  ports  are  fitted  with  Individually  fitted 
plugs,  which  can  be  removed  for  insertion  of  probe  traverse  mechanisms, 
etc.  Several  static  pressure  taps  are  provided  on  the  top  surface  of 
the  channel  throughout.  Both  the  probe  ports  and  the  pressure  taps  are 
spaced  7.6  cm  along  the  channel  length. 

The  development  section,  the  test  section,  the  recovery  section, 
and  the  exit  section  are  constructed  of  1.3  an  thick  plexiglas  rein¬ 
forced  at  periodic  intervals  by  plexiglas  ribs.  The  use  of  transparent 
plexiglas  permitted  easy  optical  access  for  the  application  of  laser 
anemometry. 

In  order  to  permit  LUA  measurements  close  to  the  wall,  the  top  wall 
around  the  measurement  location  In  the  test  section  is  deflected 
slightly  Into  the  channel.  A  large  bolt,  extending  from  a  rigid  alumi¬ 
num  frame  placed  around  the  measurement  location,  is  used  to  press  down 
on  the  top  wall  and  deflect  it  inward. 


Significant  errors  in  the  near-wall  unsteady  measurements  can 
result  from  periodic  variation  at  the  exciting  frequency  of  the  probe 
location  with  respect  to  the  wall  (see  Acharya  and  Reynolds  L1975J, 
Simpson  et  al.  L19B3J).  The  bolt  and  frame  combination  mentioned  above 


served  the  additional  purpose  of  stiffening  the  top  Mall  of  the  channel 
around  the  measurement  location  and  therefore  reducing  any  oscillations 
of  the  top  Mall  during  unsteady  operation.  From  a  careful  evaluation  of 
the  near-Mali  unsteacty  data  (see  Appendix  B),  It  Is  believed  that  the 
present  experiments  do  not  have  significant  errors  due  to  relative  move¬ 
ment  of  the  LDA  measuring  volume  Mlth  respect  to  the  top  Mall. 

3.1.3  Oscillation  Mechanism 

Water  that  enters  the  test  section  exits  the  system  through  tMO 
different  paths.  Some  of  It  Is  bled  off  through  the  porous  bottom  Mall 
of  the  test  section.  The  rest  of  It  passes  through  the  recovery  section 
and  exits  via  the  porous  plate  In  the  exit  section.  These  paths  lead  to 
tMO  adjacent  but  separated  ducts.  A  plate  Mlth  a  number  of  slots  acts 
as  a  gate  valve  that  controls  the  exit  of  both  ducts.  There  Is  a  spec¬ 
ially  fabricated  Teflon  gasket  betMeen  the  exit  of  the  ducts  and  one  of 
the  sides  of  the  gate  valve.  The  gasket  serves  as  a  Mater  seal.  The 
other  side  of  the  gate  valve  Is  open  to  the  atmosphere.  Water  that 
passes  through  the  slots  In  the  gate  valve  drains  doMn  Into  the  sump. 

The  gate  valve  can  slide  back  and  forth  freely  on  ball-bearing 
supports.  The  slot  area  presented  to  each  duct  varies  as  the  valve 
moves,  but  the  total  slot  area  seen  by  both  ducts  is  alMays  a  constant. 
This  area  Is  the  main  controlling  resistance  for  the  floM  of  Mater  from 
the  overhead  tank  to  the  sump  beloM.  Since  It  Is  a  constant,  the  total 
floM  through  the  valve  Is  a  constant  Independent  of  the  valve  position. 
Thus  the  floM  rate  through  the  system  Is  independent  of  the  position  of 
the  valve. 

The  fraction  of  the  total  slot  area  seen  by  the  duct  connected  to 
the  test  section  Is  a  linear  function  of  the  valve  position.  Thus,  to 
the  first  order,  the  proportion  of  the  total  floM  bled  off  In  the  test 
section  Is  a  linear  function  of  the  valve  position.  The  gate  valve  Is 
oscillated  back  and  forth  sinusoidally  by  a  Skotch-yoke  mechanism  driven 
by  a  constant-speed  DC  motor.  This  causes  the  fraction  of  the  total 
f Iom  bled  off  In  the  test  section  to  oscillate  nearly  sinusoidally. 
Since  the  suction  through  the  porous  plate  Is  distributed  throughout  the 
length  of  the  test  section,  the  tree-stream  velocity  decreases  nearly 
linearly  Mlth  X  and  oscillates  In  response  to  the  oscillating  valve. 

33 


To  compensate  for  the  required  pressure  drop  across  the  bottom-wal 1 
porous  plate  in  the  test  section,  extra  corrections  slots  that  are  ex¬ 
posed  only  during  some  portions  of  the  oscillation  cycle  are  used  in  the 
gate  valve.  These  slots  are  carefully  trimmed  for  each  of  the  oscilla¬ 
tion  amplitudes,  to  minimize  the  residual  variation  in  the  flow  rate 
upstream,  as  the  valve  oscillates. 

The  system  is  carefully  adjusted  so  that  the  free-stream  velocity 
remains  a  constant  throughout  the  test  section  when  the  valve  is  kept 
stationary  at  the  position  corresponding  to  the  0°  cycle-phase  angle. 
By  varying  the  speed  of  the  DC  motor,  the  frequency  of  test-section 
free-stream  velocity  oscillations  can  be  varied  over  a  range  of  0  to  2 
Hz.  By  varying  the  amplitude  of  the  gate  valve  oscillations,  the 
amplitude  of  free-stream  velocity  oscillations  can  be  varied.  The 
maximum  usable  free-stream  velocity  amplitude  at  the  end  of  the  test 
section,  normalized  by  the  mean  free  stream  velocity  at  the  entrance  is 
about  30%.  This  limit  is  set  by  the  excessive  growth  of  the  boundary 
layer  at  low  frequencies,  for  large  free-stream  adverse  gradients,  in 
comparison  to  the  channel  height. 

3.1.4  Water  Treatment  and  Control 

Tap  water  contains  a  large  number  of  microscopic  particles.  For 
proper  operation  of  the  Laser  Anemometer  system,  it  was  found  necessary 
to  filter  the  water  to  remove  particles  larger  than  about  5  microns.  In 
order  to  maintain  a  constant  flow  temperature  over  long  periods,  cooling 
is  continuously  provided  by  passing  part  of  the  flow  through  a  water- 
chiller  unit.  A  thermistor  located  in  the  exit  section  is  used  to  sense 
the  water  temperature,  and  a  closed- loop  electronic  controller  is  used 
to  vary  the  duty  cycle  of  the  chiller.  With  this  arrangement,  the  water 
temperature  can  be  maintained  to  within  0.1°C.  All  experiments  were 
conducted  at  a  nominal  temperature  of  17°C. 

Chemical  treatment  of  the  water  to  prevent  algae  growth  was  not 
possible,  since  the  chemicals  available  for  this  purpose  caused  excess¬ 
ive  numbers  of  air  bubbles  to  remain  In  the  water,  rise  to  the  top  wall, 
and  prevent  near-wall  measurements.  The  relatively  low  water  tempera¬ 
ture  was  therefore  important  in  minimizing  algae  growth. 
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3 . 1  -  b  Li) A  Mount 


The  sending  and  receiving  optics  of  the  laser  anemometer  system  are 
mounted  on  two  arms  of  a  U-shaped  rigid  structure  constructed  out  of 
aluminum  channels.  These  arms  hang  on  either  side  of  the  channel.  The 
structure  is  mounted  on  a  machine  slide  with  fine  traverse  capability  in 
the  Y  (vertical)  direction.  The  slide  is  attached  to  a  carriage  that 
can  move  in  the  L  (spanwise)  direction  on  precision  guide  rods.  The 
guide  rods  are  mounted  on  a  platform  that  can  move  in  the  X  direction 
on  wheels  supported  by  I-beams  that  run  parallel  to  the  channel.  The 
platform  can  be  locked  in  a  given  place  to  the  I-beams.  Since  the  col¬ 
lecting  and  receiving  optics  are  fixed  to  the  same  rigid  structure,  the 
alignment  of  the  optics  is  not  affected  as  the  system  is  moved  about. 


Data-Acquisition  and  Control  System 


Much  of  the  experimental  data  acquisition  and  control  procedure  was 
automated  to  meet  the  present  research  objectives  in  a  reasonable  amount 
of  time.  The  typical  data-acquisition  times  were  very  long.  For  a  good 
profile  definition,  a  large  number  of  points  (3U-5U)  was  needed  in  the 
boundary  layer.  For  a  good  signal-to-noise  ratio  in  the  phase-averaged 
velocity  and  turbulence  (especially  in  the  low-amplitude  experiment),  a 
large  number  of  cycles  (500)  was  needed.  For  0.1  Hz,  the  above  require¬ 
ments  translate  to  over  100  hours  of  data  collection.  The  frequency  of 
oscillation  needed  to  be  maintained  constant  over  such  long  times. 
Several  frequencies  and  X  stations  needed  to  be  studied  to  provide  a 
good  understanding  of  the  flow  behavior  and  to  provide  an  extensive  data 
set.  A  computer-based  and  automated  data-acquisition  and  control  system 
was  ideally  suited  to  meet  the  above  requirements. 


The  data-acquisition  and  control  system  (Fig.  3.3)  consists  of  the 
laboratory  computer  system,  the  Laser-Doppler  Anemometer  system,  the 
shaft  position  encoder  system,  the  LjA  traverse  system,  the  frequency 
control  system,  and  the  pressure-measurement  system.  Each  of  these  sub¬ 
systems  will  be  briefly  described  below. 
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3.2.1  Laboratory  Computer  System 


The  data-acqulsltlon  and  control  system  is  based  on  a  MINC  micro¬ 
computer  system  (Fig.  3.4)  consisting  of  a  Digital  Equipment  Corporation 
LSI  11/2  processor  with  64  KB  of  memory.  Two  single-sided  double¬ 
density  diskette  drives  provide  1  MB  of  mass  storage.  The  system  is 
equipped  with  a  VT105  Graphics  Terminal  for  interactive  usage,  an  IUS460 
Paper  Tiger  Matrix  Printer,  and  a  HP9872C  Graphics  Plotter  interfaced 
through  an  IEEE-488  Instrument  Bus  controller.  The  system  has  been 
configured  with  the  following  hardware  to  facilitate  interfacing  witn 
laboratory  instruments  and  control  systems  for  real-time  data  acquisi¬ 
tion  and  control: 

Real  Time  Clock:  The  clock  can  be  programmed  to  provide  periodic 
interrupts  from  0.002  Hz  to  10^  Hz  and  to  time  external  events;  the 
clock  also  permits  synchronization  with  external  events. 

Analog-to-Digital  Converter:  The  A/D  has  eight  fully  differential 
channels,  a  programmable  gain  amplifier  (gains  of  1,  2,  4,  and  8), 
12-bit  accuracy,  35  KHz  conversion  rate,  and  a  swich-selectable  range 
of  (0  to  10)  V  or  (-10  to  +10)  V. 

Digital-to-Analog  Converter:  The  D/A  has  four  channels  of  12-bit 
accuracy  and  switch- selectable  ranges  of  (0  to  5)  V,  (0  to  10)  V, 
or  (-10  to  +10)  V. 

Digital  Input  and  Output:  Two  channels  of  DIN  and  two  channels  of 
DOUT  are  provided.  Each  channel  Is  16-bits  wide  and  Includes  lines  for 
triggering  and  sensing  service  requests  from  external  devices. 

All  the  Input  and  output  connectors  and  switches  of  the  above  in¬ 
terface  boards  are  brought  to  the  front  panel  of  the  computer  for  easy 
connection  to  laboratory  Inst.iments. 

3.2.2  Laser  Doppler  Anemometer  System 

The  LDA  system  Is  a  one-component  forward- scatter,  differential, 
Doppler-mode  system  employing  frequency  shifting.  The  Doppler  signal  is 
processed  by  a  frequency  tracker  working  on  the  phase- locked  loop  prin¬ 
ciple.  The  system  is  composed  of  the  following  DISA  55X  Modular  LDA 
Optical  components. 


5  mW  He-Ne  coaxial  laser  head, 

55X24  beam  splitter, 

55X29  Bragg  cell  (40  MHz  frequency  shift), 

55X28  beam  displacer, 

55X52  front  lens  (focal  length  600  mm), 

55X34  photo-multiplier  optics, 

55X08  photo-multiplier  section, 

55N10  frequency  shifter, 

55N20  frequency  tracker. 

The  sending  and  receiving  optics  are  mounted  on  two  rigid  arms  of  the 
laser  mount  structure  that  extend  down  on  either  side  of  the  water 
channel.  Two  optical-quality  mirrors  are  used  for  turning  the  beams 
into  the  water  channel  and  for  focusing  of  the  collecting  optics.  The 
frequency  shifter  unit  can  electronically  vary  the  frequency  (around  4U 
MHz)  by  which  one  of  the  beams  is  shifted  with  respect  to  the  other. 
The  frequency  tracker  has  both  a  front-panel  digital  readout  of  the 
velocity  and  an  analog  voltage  output  proportional  to  the  instantaneous 
tracking  frequency,  which  in  turn  is  proportional  to  the  average  veloc¬ 
ity  of  particles  in  the  measuring  volume  at  that  instant. 


3.2.3  Shaft  Encoder  System 

The  shaft  encoder  system  consists  of  a  Litton  Model  76  Absolute 
Shaft  Position  Encoder  and  the  associated  electronics.  The  encoder  is 
connected  to  the  shaft  of  the  UC  motor  that  drives  the  oscillation  mech¬ 
anism  and  provides  a  10-bit  Gray  coded  output  of  the  angular  position  of 
the  shaft.  The  0°  output  of  the  encoder  is  made  to  correspond  to  the  u° 
cycle-phase  angle.  The  encoder  electronics  provides  the  power  needed  to 
operate  the  encoder  and  converts  the  Gray  code  to  natural  binary  code 
for  input  to  the  computer.  In  addition,  the  encoder  electronics  pro¬ 
vides  sampling  triggers  at  512  equally  spaced  phases  in  the  cycle  and  a 
pulse  at  the  occurrence  of  the  0°  cycle-phase  angle. 


3.2.4  LDA Mount  Traverse  System 


A  SLO-SYN  Stepping  Motor  (Model  M093-FC07)  is  geared  to  the  shaft 
of  the  slide  of  the  LDA  Mount  system.  The  motor  Is  energized  by  a 
SLO-SYN  Translater  (Type  ST101).  The  translater  Is  controlled  by  an 
electronic  Probe  Controller.  The  probe-controller  electronics  can  be 
operated  from  the  front  panel  or  can  be  controlled  by  the  computer.  The 
controller  can  be  set  to  provide  digital  LED  readouts  of  the  Y  loca¬ 
tion  with  respect  to  the  wall  of  the  measuring  volume  of  the  LDA  In 
units  of  0.0025  cm.  Each  step  of  the  stepping  motor  corresponds  to 
0.00032  cm  movement  of  the  measuring  volume. 

3.2.5  Pressure  Measurement  System 

The  pressure  measurement  system  Is  based  on  a  48-channel  scanlvalve 
pressure  switch  and  a  Vail  dyne  Model  MP245  pressure  transducer.  Pres¬ 
sure  transmitting  lines  from  all  the  pressure  taps  In  the  water  channel 
are  brought  to  the  Scanlvalve.  There  Is  provision  for  purging  these 
lines  of  trapped  air  bubbles.  The  pressure  switch  is  capable  of  rapid 
scanning  of  the  channels  without  Imposing  a  volume  change  on  the  trans¬ 
ducer.  The  transducer  cavities  can  also  be  purged  of  air  bubbles. 

The  range  of  differential  pressure  that  the  transducer  can  measure 
is  from  0  to  10  cm  H20.  The  actual  range  over  which  the  transducer  was 
used  In  the  present  experiments  was  from  2  to  3  cm  H20.  Over  this 
range,  the  transducer  and  the  associated  electronics  were  calibrated  to 
an  accuracy  better  than  0.1%.  Two  columns  of  water  whose  height  differ¬ 
ence  could  be  measured  to  better  than  0.025  mm  were  used  to  provide  the 
calibration  differential  pressure  for  the  transducer. 

The  pressure  measurement  system  was  used  In  the  present  work  only 
for  measuring  the  free-stream  velocity  In  the  channel  under  steady 
operating  conditions,  and  for  the  calibration  of  the  LDA  system  (see 
Section  3.3.3).  The  free-stream  velocity  could  be  measured  to  an  accu¬ 
racy  of  better  than  0.1%. 


3.3  Data-Acquisition  and  Control  Procedures 

The  data-acqulsition  and  control  procedures  deal  with  the  LUA  probe 
volume  positioning,  the  calibration  of  the  LDA  system,  the  long-term 
control  of  frequency,  and  the  experimental  run  procedure.  Each  of  these 
is  described  below. 

3.3.1  LUA  Probe  Volume  Positioning 

The  LDA  traverse  system  permits  very  precise  movement  (0.02b  an;  of 
the  LDA  probe  volume  in  the  Y  direction.  However,  each  time  the  LDA 
system  is  moved  to  a  different  X  location,  it  is  necessary  to  position 
the  measuring  volume  accurately  with  respect  to  the  wall.  The  means  and 
the  procedure  for  doing  this  are  described  below. 

A  standard  tunnel  port  plug  was  modified  to  accommodate  at  its  cen¬ 
ter  a  stainless  steel  stem  with  a  disk  on  the  end  toward  the  inside  of 
the  tunnel.  The  stem  is  spring-loaded  so  that,  in  its  normal  retracted 
position,  the  disk  would  be  flush  with  the  wall  of  the  tunnel.  The  stem 
and  the  disk  can  be  pushed  into  the  tunnel  and  locked  in  place.  A  slot 
is  cut  on  the  stem  surface  perpendicular  to  the  stem  axis.  A  sharp 
needle  is  cemented  in  the  slot.  The  distance  of  the  needle  from  the 
surface  of  the  plug  is  carefully  measured  using  an  optical  comparator. 

The  LDA-mount  carriage  is  moved  to  the  desired  X  location  and 
locked  in  place.  The  LDA-mount  slide  is  moved  so  that  the  probe  volume 
intersects  the  tip  of  the  needle.  The  needle  casts  a  sharp  shadow  be¬ 
hind,  and  the  bright  streaks  of  laser  light  on  either  side  of  the  shadow 
are  made  equal  in  extent  in  order  to  position  the  center  of  the  probe 
volume  with  the  tip  of  the  needle.  This  positioning  is  accurately 
repeatable  to  within  0.02b  mm. 

3.3.2  Frequency  Control 

The  DC  motor  controller  is  capable  of  keeping  the  frequency  of 
oscillation  constant  to  within  0.2%  in  the  short  term.  Over  the  long 
term,  however,  due  to  thermal  drift  of  analog  components  used  in  the 
controller,  the  frequency  can  vary  considerably.  This  is  especially 
acute  for  low  frequencies,  since  the  total  data-acqulsition  time  needed 
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for  a  given  number  of  cycles  Is  large.  The  problem  was  overcome  by 
controlling  the  set-point  voltage  for  the  controller  externally  by  the 
computer. 

The  transfer  function  from  the  external  DC  voltage  Input  of  the 
motor-speed  controller  to  the  oscillation  frequency  was  found  to  be  very 
highly  linear.  One  channel  of  the  D/A  Is  connected  to  this  input. 
First,  the  D/A  output  is  set  at  the  nominal  voltage  corresponding  to  the 
desired  frequency,  and  the  real-time  clock  is  set  up  to  time  external 
events.  The  clock  is  started  on  an  occurrence  of  the  0°  position  and 
stopped  on  the  next  occurrence.  The  output  of  the  clock  is  then  the 
period  of  oscillation  from  which  the  actual  frequency  is  calculated. 
The  error  in  the  actual  frequency  is  used  to  correct  the  D/A  output 
voltage  appropriately  and  the  procedure  repeated  until  the  actual  and 
the  desired  frequencies  differ  by  less  than  0.2%.  With  this  procedure, 
it  was  possible  to  keep  the  oscillation  frequency  constant  to  within 
0.2%  over  any  desired  length  of  time. 

3.3.3  LDA  Calibration 

An  LDA  system  is  a  linear  transducer  of  velocity  to  Doppler  fre¬ 
quency.  Ideally,  if  the  separation  of  the  beams,  the  front  lens  focal 
length  and  the  details  and  properties  of  the  media  along  the  beam  path 
are  known,  then  the  proportionality  constant  can  be  easily  calculated. 
In  practice,  lack  of  precise  knowledge  of  any  of  the  parameters  makes  it 
desirable  to  do  a  one-point  calibration  using  another  means  of  indepen¬ 
dent  measurement  of  velocity.  Certainly,  it  is  advisable  to  cross-check 
any  computed  calibration  factor  with  a  measured  one.  The  procedure 
adopted  in  the  present  case  is  described  below. 

The  LDA  measuring  volume  and  a  pitot  tube  were  positioned  in  the 
free  stream.  The  pressure  difference  between  the  pitot  tube  and  a 
static  pressure  tap  at  the  same  X  location  was  measured  very  accu¬ 
rately  using  the  pressure  transducer,  and  the  free-stream  velocity  was 
calculated.  From  the  average  Doppler  frequency  measured  by  the  LDA  and 
the  free-stream  velocity,  the  calibration  factor  for  the  LDA  was  easily 
calculated.  This  procedure  was  repeated  over  a  20%  range  of  free-stream 
velocities,  and  the  results  were  averaged.  The  measured  calibration 
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factor  was  used  In  all  data  reduction.  An  estimate  of  the  calibration 
factor  made  from  the  beam  separation,  the  focal  length  of  the  front 
lens,  and  the  refractive  index  of  water  is  only  1.3%  different  from  the 
measured  value. 

When  the  channel  is  set  up  for  uniform  free- stream  velocity 
throughout  the  rig,  the  LDA  measurements  of  free-stream  velocity  along 
the  channel  are  within  U.5%  of  the  values  given  by  pressure  measure¬ 
ments.  This  indicates  that  the  LDA  calibration  factor  is  good  to  within 
0.5%  throughout  the  channel. 

3.3.4  Experimental  Run  Procedure 

Figure  3.3  shows  a  schematic  of  the  interconnections  among  the  dif¬ 
ferent  parts  of  the  data-acquisition  and  control  system.  Only  a  very 
brief  listing  of  the  different  steps  involved  in  running  the  experiment 
will  be  given  below.  Much  of  the  procedure  is  automated,  and  extensive 
software  was  developed  for  data  acquisition,  control,  and  processing. 
The  listing  of  all  programs  relevant  to  the  running  of  the  experiment 
can  be  found  in  Jayaraman  L1982}.  The  program  listings  should  be  con¬ 
sulted  for  the  exact  procedural  details. 

The  LDA  is  moved  to  the  desired  X  location,  and  the  probe  volume 
is  positioned  as  described  above.  The  desired  frequencies,  the  Y 
locations,  the  number  of  cycles  to  average  over,  and  the  run  title  are 
input  to  the  data-acquisition  program.  The  program  moves  the  LDA  to  a 
desired  Y  location  (starting  from  the  wall  and  moving  out),  sets  up  a 
desired  frequency,  and  starts  collecting  data.  Every  ten  cycles,  it 
verifies  the  frequency  of  oscillation  and  re- synchronizes  the  data  col¬ 
lection  with  the  0°  cycle-phase  angle.  After  data  have  been  collected 
over  the  desired  number  of  cycles,  the  program  stores  the  data  on  the 
disk,  sets  up  the  next  frequency,  and  repeats  the  process  until  all 
frequencies  are  done  at  that  Y  location.  The  program  then  moves  the 
LDA  to  the  next  Y  location  and  repeats  the  process  until  all  Y 
locations  are  done. 

The  whole  system  has  been  reliably  operated  In  such  an  automated 
mode  over  periods  as  long  as  240  hours  with  very  minimal  operator  inter- 
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Extensive  software  was  developed  for  the  purposes  of  analyzing  the 
data  obtained  by  the  data-acqulsltlon  programs.  Software  packages  for 
plotting  of  the  results  using  the  plotter  were  also  developed.  All 
these  programs  are  listed  In  Jayaraman  L1982J  and  should  be  consulted 
for  precise  details  of  the  data-reductlon  and  plotting  procedures. 

The  data-acqulsltlon  programs  obtained  phase- averaged  streamwise 
velocity  component  (<  u  >)  and  its  square  (<  u^  >).  From  these  the 
mean  and  the  periodic  components  of  the  streamwise  velocity  component 
and  the  square  of  streamwise  turbulent  fluctuations  were  calculated  as 
shown  below. 


U  «  <  u  >  ;  u  -  <  u  >  -  U  ; 

<  u'u'  >  ■  <  u*  >  -  <  u  >*  ; 


(u'u‘)  ■  <  u'u'  >  ;  (u’u')  *  <  u'u'  >  -  u'u' 

The  mean  velocity  data  were  analyzed  using  the  procedure  described 
by  Coles  [1968]  to  determine  the  wall  shear  stress  and  the  boundary 
layer  thicknesses.  The  boundary  layer  thickness  so  determined  will  be 
referred  to  as  the  "Coles  boundary  layer  thickness"  In  the  rest  of  the 
report.  It  Is  a  much  more  definite  parameter  than  the  99%  thickness. 

The  periodic  component  of  the  velocity  and  the  turbulence  were 
analyzed  using  an  FFT  procedure  to  determine  the  amplitude  and  the  phase 
of  the  first  harmonic.  The  variation  over  the  cycle  of  displacement  and 
momentum  thickness  and  shape  factor  was  computed  easily  from  the  phase- 
averaged  velocity  profiles. 
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Chapter  4 

EXPERIMENTAL  RESULTS 


Two  major  series  of  experiments  were  done.  In  the  first,  the 
amplitude  of  free- stream  oscillations  at  the  end  of  the  test  section 
normalized  by  the  upstream  free-stream  mean  velocity  was  about  5%;  in 
the  second,  it  was  about  25%.  These  experiments  will  be  referred  to  as 
the  "low-amplitude"  and  the  "high-amplitude"  experiment,  respectively. 
In  both  experiments,  data  were  collected  at  several  different  frequen¬ 
cies  and  X  stations.  Tables  4.1  and  4.2  summarize  all  the  cases 
covered  by  the  two  experiments  and  Include  the  corresponding  Strouhal 
numbers,  defined  as  follows: 


4X-XJ 


c  0)6 

s.  *  — 

6  U 


where 


X  is  the  streamwlse  distance  at  the  measurement  station, 

XQ  is  the  streamwise  distance  of  the  entrance  to  the  test  section, 

u  is  the  circular  frequency  in  radians/ second, 

Ug  Is  the  local  free-stream  mean  velocity,  and 

6  is  the  thickness  of  the  mean  boundary  layer  determined  by  Coles 
analysis  procedure  (Section  3.4). 

The  streamwlse  distance  in  the  test  section  is  normalized  by  the  test- 
section  length  as  follows: 

(X-X  ) 

u  •  w 


where  L  Is  the  length  of  the  test  section,  60.9  cm.  It  can  be  seen 
that  the  experiments  cover  a  very  wide  range  In  Strouhal  number. 

Detailed  results  for  all  the  cases  studied  are  available  In  a  sep¬ 
arate  report  [Jayaraman,  Parikh,  and  Reynolds,  1982]  on  microfiches. 
Refer  to  Appendix  A  for  the  organization  and  contents  of  the  micro¬ 
fiches.  In  the  following  presentation  of  the  results  of  the  experi¬ 
ments,  only  representative  cases  will  be  used.  They  will  be  primarily 
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from  the  high-amplitude  experiment  supplemented  by  the  low-amplitude 
cases,  where  needed  to  show  the  effect  of  amplitude.  If  any.  The  term 
"turbulence"  Is  used  to  denote  the  square  of  streamwlse  turbulent  fluc¬ 
tuations  (u'u1). 

4.1  General  Considerations 

In  evaluating  the  results  to  be  presented,  several  Important  con¬ 
siderations  should  be  kept  In  mind.  These  are  addressed  below. 

4.1.1  Number  of  Cycles 

In  all  cases,  phase-averaged  streamwlse  velocity  and  the  square  of 
streamwlse  turbulent  velocity  fluctuations  were  computed  over  500  cycles 
of  oscillation.  The  number  of  cycles  chosen  was  a  compromise  between 
very  long  run  times  and  good  slgnal-to-nolse  ratios  In  the  averages. 
Naturally,  the  high-amplitude  cases  have  much  better  slgnal-to-nolse 
ratios  than  do  the  low-amplitude  cases.  For  this  reason,  the  hlgh- 
amplltude  experiment  Is  more  extensive  than  the  low-amplitude  one. 

4.1.2  Zero  Frequency  Measurements 

There  Is  one  Important  difference  between  the  two  experiments.  In 
the  case  of  the  low-amplitude  experiment,  the  flow-oscillation  gate 
valve  could  be  kept  stationary  at  different  locations,  and  the  resulting 
steady  velocity  profiles  could  be  measured  and  used  for  constructing  the 
response  of  the  boundary  layer  when  subjected  to  extremely  low  frequen¬ 
cies  of  oscillation.  This  could  not  be  done  In  the  high-amplitude 
experiment,  since  under  large  adverse  free-stream  gradient  conditions 
the  boundary  layer  separated  catastrophically  from  the  wall  and  a  vis¬ 
cous  region  filled  up  the  entire  channel.  However,  It  was  found  that 
flow  oscillation  at  frequencies  as  low  as  0.1  Hz  was  sufficient  to  pre¬ 
vent  the  boundary  layer  from  such  separation,  and  meaningful  unsteady 
data  could  be  collected. 
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4.1.3  Very  Near-Wall  Data  Restrictions 

The  size  of  the  LDA  probe  volume  In  the  Y  direction  Is  estimated 
to  be  approximately  0.3  mm.  The  size  of  the  actual  measuring  volume  is 
determined  by  the  probe  volume  and  the  collecting  optics,  and  is  likely 
to  be  of  the  order  of  0.3  mm.  The  output  signal  of  the  photomultiplier 
is  an  average  of  signals  from  particles  in  different  parts  of  the  mea¬ 
suring  volume.  If  the  probe  volume  center  is  moved  closer  to  the  wall 
than  its  radius,  then  part  of  the  volume  is  cut  off  by  the  wall.  The 
effective  center  of  such  a  truncated  volume  Is  farther  from  the  wall 
than  the  probe  volume  center.  It  is  estimated  that  data  closer  to  the 
wall  than  0.2  mm  suffer  from  the  above  restrictions,  due  to  finite 
measuring  volume  size. 

The  effect  on  mean  velocity  is  to  measure  higher  velocity  than 
actual  at  a  given  distance  from  the  wall.  Another  way  to  look  at  this 
is  to  say  that  the  measured  velocity  is  to  be  ascribed  to  a  location 
from  the  wall  farther  than  the  nominal  location.  Data  from  such  loca¬ 
tions  are  still  included  in  the  results  to  be  presented  below.  It  is 
believed  that  they  are  useful  for  extracting  qualitative  trends.  Ho 
conclusions  will  be  drawn  based  solely  on  data  from  such  locations 
without  attendant  warnings  to  that  effect. 

It  should  be  mentioned  that  the  effect  of  these  points  near  the 
wall  on  computation  of  boundary  layer  integral  thicknesses  and  shape 
factor  is  very  small.  This  can  be  seen  by  comparing  the  results  of  com¬ 
puting  these  parameters  for  the  mean  profiles  by  two  different  methods. 
In  one  method,  standard  wall  functions  are  used  up  to  a  Y*  of  50  and 
the  experimental  data  beyond,  as  suggested  by  Coles  [1968].  In  the 
other  scheme,  only  experimental  data  are  used  throughout.  The  differ¬ 
ence  between  the  two  methods  in  all  cases  is  found  to  be  less  than  0.2%. 

4.1.4  Estimated  Uncertainties 

Recall  that  one  of  the  main  objectives  of  this  work  is  to  provide  a 
definitive  data  set  on  a  well-defined  unsteady  turbulent  boundary  layer 
for  the  purposes  of  evaluating  the  performance  of  computational  methods 
and  turbulence  models  for  such  flows.  As  pointed  out  by  Moffat  [1981], 
such  a  comparison  is  meaningful  only  in  the  light  of  the  Nth-order 
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uncertainties  in  the  measured  and  computed  results.  The  Ntn-order 
uncertainty  in  a  measured  value  is  an  estimate  of  Its  absolute  accuracy 
and  includes  the  scatter  on  repeated  trials  (lst-order  uncertainty) . 

The  estimated  N^-order  uncertainties  in  the  measured  and  calcula¬ 
ted  quantities  are  shown  in  Table  4.3.  The  associated  confidence  level 
Is  20  to  1.  The  Nth-order  uncertainties  were  propagated  using  the  pro¬ 
cedure  set  forth  by  Kline  and  McClintock  [1953].  Note  that  many  of  the 
estimates  are  dependent  on  the  location  in  the  boundary  layer  and  the 
amplitude  of  oscillation. 

4.2  Flow  Qualification 

The  purpose  of  this  section  is  to  establish  that  the  flow  in  the 
test  rig  is  spanwise  uniform,  the  growth  of  the  boundary  layer  in  the 
rig  under  zero  pressure-gradient  conditions  is  normal,  and  the  induced 
disturbance  upstream  during  unsteady  operation  is  acceptably  small. 
Each  of  these  issues  is  discussed  below. 


4.2.1  Spanwise  Uniformity 

To  determine  the  two-dimensionality  of  the  flow  in  the  tunnel, 
closely  sampled  spanwise  surveys  of  velocity  at  a  fixed  location  from 
the  wall  in  the  boundary  layer  were  taken  by  moving  the  LDA  measuring 
volume  in  the  spanwise  direction.  In  the  pre-test  operation  of  the  test 
rig,  with  60%  open  area  ratio  screens,  it  was  found  that  the  mean  veloc¬ 
ity  varied  periodically  in  the  spanwise  direction,  with  the  peak-to-peak 
variation  being  more  than  10%  of  the  free- stream  velocity.  The  cause  of 
this  variation  was  ultimately  traced  to  the  screens  upstream.  It  ap¬ 
pears  that  the  screen  surface  should  be  stretched  tight  to  eliminate  any 
wrinkles  that  cause  longitudinal  vorticity,  which  causes  spanwise  vari¬ 
ation  in  velocity.  Reconstruction  of  the  screen  box  with  plexiglas 
permitted  very  easy  detection  of  deterioration  of  cleanliness  of  the 
screens.  Use  of  70%  open  area  ratio  screens  that  can  be  removed  period¬ 
ically  for  cleaning  and  kept  very  tightly  stretched  in  place  put  the 
problem  under  complete  control. 


Figures  4-la-d  show  typical  spanwise  variations  In  the  mean  veloc¬ 
ity  (normalized  by  the  mean  free-stream  velocity  at  the  entrance  to  the 
test  section,  UQ) ,  the  mean  turbulence  (normalized  by  the  square  of 
UQ),  the  amplitude  of  velocity  oscillations  (normalized  by  UQ)  and 
its  phase  angle  at  X'  *  U.94  for  the  low-amplitude  experiment.  Fig¬ 
ures  4.1e-f  show  the  same  for  the  high-ampl itude  experiment.  The  span- 
wise  uniformity  is  worst  at  the  end  of  the  test  section,  and  this  point 
was  chosen  for  that  reason.  The  abscissas  are  normalized  by  the  channel 
width.  Note  that  the  data  points  cover  the  central  half  of  the  channel 
and  that  they  are  closely  spaced  (about  a  quarter  of  the  mean  boundary 
layer  thickness). 

The  mean  velocity  is  constant  to  within  1%  across  the  span,  while 
the  root-mean-squared  turbulence  is  uniform  to  within  1.5%.  The  vari¬ 
ation  in  the  amplitude  of  velocity  oscillations  is  within  0.3%,  while 
that  of  the  phase  angle  is  less  than  3°  across  the  span.  The  worst 
spanwise  nonuniformity  is  found  in  the  mean  turbulence  in  the  high- 
amplitude  case.  It  is  interesting  that  this  variation  does  not  show  up 
in  the  mean  velocity,  the  amplitude,  or  the  phase. 

The  results  presented  are  at  one  fixed  location  from  the  wall  for 
one  frequency  of  oscillation.  Measurements  at  other  Y  locations,  X 
stations,  and  frequencies  show  variation  of  the  same  order.  Thus  it  can 
be  concluded  that  the  flow  in  the  central  half  of  the  channel  is  very 
highly  two-dimensional  with  respect  to  both  the  mean  and  the  periodic 
flow  fields. 


4.2.2  Zero  Pressure-Gradient  Boundary  Layer  Development 


Figure  4.2  shows  the  free-stream  velocity  variation  with  respect  to 
X  in  the  test  rig  when  the  flow-oscillation  gate  valve  is  set  at  the 
position  corresponding  to  0°  cycle-phase  angle.  Note  that,  with  careful 
adjustment  of  the  bottom-wall  suction  in  the  development  section  to 
account  for  the  blockage  due  to  the  growth  of  boundary  layers,  the  free- 
stream  velocity  is  uniform  to  within  0.3%.  This  adjustment  needed  very 
little  changes  over  very  long  periods  of  operation  of  the  rig. 


Figure  4.3  shows  the  development  of  the  mean  velocity  in  the  devel¬ 
opment  section  in  inner  coordinates.  For  comparison,  the  corresponding 


Coles  velocity  profiles  are  also  shown  In  the  figure.  Note  the  normal 
log-linear  region  and  the  growth  of  the  wake  component  as  the  boundary 
layer  grows. 

Figures  4.4a-b  show  the  development  of  the  Coles  boundary  layer 
thickness  and  the  displacement  thickness  under  zero  pressure  gradient 
conditions  In  the  rig.  Both  grow  as  a  power  of  X  somewhat  less  than  1 
(-  0.8). 

Figure  4.5  shows  the  skin-friction  coefficient  variation  with  re¬ 
spect  to  momentum  thickness  Reynolds  number  compared  to  results  from 
Anderson  et  al.  [1972]  for  the  same  range  of  Reynolds  numbers.  The 
agreement  Is  quite  good. 

Figure  4.6  shows  the  defect  velocity  profile  at  the  last  measurable 
X  station  In  the  upstream  development  section.  The  data  are  compared 
with  data  from  Klebanoff  [1955],  and  the  agreement  is  excellent.  The 
mean  turbulence  profile  In  the  same  boundary  layer  Is  compared  to  Kleba- 
noff's  data  In  Fig.  4.7.  The  agreement  is  very  good  in  the  outer  part 
of  the  boundary  layer,  but  the  present  data  are  somewhat  lower  in  the 
near-wall  region.  The  boundary  layer  momentum  thickness  Reynolds  number 
for  the  present  boundary  layer  Is  about  2500,  whereas  that  of  Klebanoff 
Is  7300.  This  Is  thought  to  be  the  cause  of  this  difference. 

A  global  check  on  the  two-dimensionality  of  the  flow  Is  to  perform 
a  momentum  balance  on  the  data.  The  procedure  followed  is  the  same  as 
that  described  In  Coles  and  Hirst  [1968].  As  pointed  out  by  Moffat 
[1981],  this  balance  is  a  useless  exercise  without  an  assessment  of  the 
significance  of  any  Imbalance,  In  light  of  the  estimated  Nth-order  un¬ 
certainty  In  closing  the  momentum  Integral  equation.  Figure  4.7.1  is  a 
plot  of  the  left-  and  the  right-hand  sides  of  the  integrated  form  of  the 
momentum  Integral  equation,  normalized  by  the  momentum  deficit  at  the 
first  profile.  Shown  also  Is  the  estimated  Nth-order  uncertainty  in 
closing  the  equation.  It  Is  clear  that  the  momentum  balance  is  within 
the  estimated  uncertainty. 

Thus  one  can  conclude  that  the  growth  of  the  boundary  layer  in  the 
test  r4g  under  zero  pressure-gradient  conditions  Is  normal  and  two- 
dimensional.  For  reference  purposes,  some  of  the  parameters  associated 


with  the  boundary  layer  at  the  entrance  to  the  test  section  are  given 
below. 


Free- stream  velocity,  U0 
Coles  boundary  layer  thickness,  6 
Displacement  thickness, 

Momentum  thickness, 

Shape  factor,  H 

Momentum  thickness  Reynolds  Number, 
Skin  friction  coefficient,  Cj 
Kinematic  viscosity,  v 


732  mm/s 
35.5  mm 
5.9  mm 
4.2  mm 
1.42 
2790 

3.22  x  10'3 
1.09  x  10-6  m2/s 


4.2.3  Upstream  Disturbance 

Ideally,  there  would  be  no  variation  in  the  upstream  flow  as  the 
flow  is  oscillated  in  the  test  section.  This  would  be  the  case  if  there 
were  no  flow  capacitance  in  the  system  and  if  the  only  resistance  to  the 
flow  were  due  to  the  slots  in  the  gate  valve.  Neither  of  the  above 
requirements  is  met  exactly  in  reality.  There  is  a  finite  flow  capaci¬ 
tance  in  the  system,  due  perhaps  to  the  flexibility  of  plexiglas.  To 
have  uniform  suction  in  the  bottom-wall  porous  plat^  in  the  test  sec¬ 
tion,  there  needs  to  be  sufficient  pressure  drop  across  the  plate,  com¬ 
pared  to  the  static  pressure  change  in  the  test  section.  Nonuniformity 
of  suction  in  the  test  section  would  result  in  phase  variation  in  the 
free  stream,  causing  the  phase-averaged  free  stream  to  be  nonlinear. 
Extra  correction  slots  that  are  exposed  only  during  some  portion  of  the 
oscillation  cycle  were  needed  in  the  gate  valve  to  compensate  for  the 
required  pressure  drop  across  the  porous  plate.  Both  the  correction 
slot  and  the  sizing  of  the  holes  in  the  porous  plate  needed  to  be  care¬ 
fully  trimmed  for  each  of  the  oscillation  amplitudes  studied.  Thus  it 
should  be  understood  that  the  actual  upstream  disturbance  and  the  im¬ 
posed  free-stream  excitation  in  the  experiments  to  be  presented  below 
were  the  best  achievable  in  the  test  facility  without  resorting  to  more 
active  control. 


Figure  4.8  shows  the  root-mean- square  value  of  the  velocity  oscil¬ 
lations  in  the  boundary  layer  at  the  last  measurable  X  station  in  the 
development  section,  expressed  as  a  percentage  of  the  free- stream  mean 
velocity.  The  reason  for  using  rms  values  is  that  the  organized  distur¬ 
bance  velocity  upstream  is  not  sinusoidal  for  some  frequencies.  Several 
frequencies  are  shown  for  the  high-amplitude  experiment.  It  is  clear 
that  the  disturbance  level  is  essentially  constant  within  the  boundary 
layer  and  that  the  level  increases  with  frequency.  Even  for  the  highest 
frequency,  it  is  only  about  1.6%.  The  disturbance  levels  for  the  low- 
amplitude  experiments  (not  shown)  are  even  lower. 

Figures  4.9  and  4.10  show  the  mean  velocity  and  mean  turbulence 
profiles  at  the  same  X  location  at  different  frequencies.  The  pro¬ 
files  show  no  effects  due  to  the  small  unsteadiness  in  the  free  stream. 

Thus  the  upstream  flow,  for  all  practical  purposes,  can  be  consid¬ 
ered  to  be  a  normal,  standard,  steady,  flat-plate  boundary  layer,  unaf¬ 
fected  by  the  downstream  flow  oscillations. 


4.3  Imposed  Free-Stream  Excitation  in  the  Test  Section 

The  design  objectives  for  the  imposed  free-stream  oscillations  ;n 
the  test  section  were  as  follows: 

•  linear  decrease  in  time-averaged  free-stream  velocity; 

•  linear  increase  in  the  amplitude  of  oscillations; 

•  linear  variation  of  phase-averaged  free-stream  velocity. 

These  requirements  implied  that  there  would  be  no  phase  variation  in  the 
free-stream  oscillations.  How  well  these  objectives  were  met  in  the 
actual  test  rig  is  discussed  below. 

Figures  4.11a-b  show  the  variation  of  the  free-stream  mean  velocity 
in  the  test  section  at  different  frequencies  for  the  high-  and  the  low- 
amplitude  experiments.  The  mean  velocity  in  the  upstream  development 
section  Is  used  for  normalization,  and  the  X  distance  In  the  test  sec¬ 
tion  Is  expressed  as  a  fraction  of  its  length.  The  mean  velocity  de¬ 
creases  linearly  with  distance,  as  per  the  design  objective. 


Figures  4.12a-b  show  the  variation  of  the  free- stream  amplitude  of 
oscillation  normalized  by  the  development  section  free-stream  mean 
velocity  for  the  two  experiments.  The  amplitude  grows  linearly  with 
distance,  as  per  the  design  objective.  For  the  high-amplitude  experi¬ 
ment,  the  slope  is  somewhat  lower  at  lower  frequencies.  The  large  var¬ 
iation  in  the  thickness  of  the  boundary  layer  over  the  oscillation  cycle 
at  lower  frequencies  is  believed  to  be  the  cause  of  the  lower  slope. 

Figures  4.13a-b  show  the  variation  of  the  free-stream  phase  angle 
with  respect  to  the  last  measured  station  in  the  recovery  section.  The 

phase  lag  in  the  upstream  locations  in  the  test  section  increases  with 

frequency.  Rapid  variations  in  phase  appear  where  the  oscillation 
amplitudes  are  relatively  small.  The  phase  variation  is  less  than  20° 
over  two- thirds  of  the  test  section  downstream  and  the  recovery  sec¬ 
tion.  These  phase  variations  are  believed  to  be  caused  by  the  non¬ 

uniformity  in  the  bottom-wall  suction  and  not  by  a  traveling  wave  inher¬ 
ent  to  the  flow.  They  are  expected  to  have  very  little  effect  on  the 
behavior  of  the  boundary  layer.  They  have  been  fully  documented  for 
inclusion  in  any  computational  scheme. 

Figures  4.14a-d  show  the  phase-averaged  free-stream  velocity  varia¬ 
tion  in  the  test  section  for  the  0.5  Hz  high-amplitude  case.  The  varia¬ 
tion  in  the  cycle  phase  angle  range  of  60  to  300°  is  linear,  as  per  the 
design  objective.  The  variation  around  0°  is  distinctly  nonlinear.  No 
such  problems  exist  in  the  low-amplitude  experiment,  as  shown  in  Figs. 
4.15a-d  for  the  0.5  Hz  case. 

Thus  it  is  clear  that  the  mean  and  the  amplitude  of  the  imposed 
free-stream  variation  in  the  test  section  closely  match  the  design 
requirements,  while  the  phase  does  not.  It  is  recommended  that  compu¬ 
tational  procedures  use  the  actual  free-stream  variation. 

4.4  Mean  Velocity  and  Turbulence 

Figures  4.16a-b  show  the  mean  velocity  profiles  at  one  X  station 
at  different  frequencies.  The  distance  from  the  wall  has  been  normal¬ 
ized  by  the  Coles  boundary  layer  thickness  of  the  mean  velocity  profile. 
The  mean  velocity  is  normalized  by  the  free-stream  mean  velocity.  No 
significant  variation  with  frequency  can  be  seen  in  the  low-amplitude 
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experiment.  In  the  high-amplitude  experiment}  a  small  variation  is  seen 
at  low  frequencies. 


Figures  4.17a-b  show  the  variation  of  mean  turbulence  profiles  at 
the  same  X  station.  The  results  are  the  same  as  for  the  mean  velocity 
profiles.  Once  again,  the  high- amplitude  experiment  low-frequency  cases 
show  a  small  variation  with  frequency. 

Thus  it  is  clear  that  the  mean  flow  is  essentially  unaffected  by 
unsteadiness  in  agreement  with  all  the  previous  works  reviewed  in  Sec¬ 
tion  1.2.  Results  from  other  X  stations  are  presented  in  Jayaraman, 
Parikh,  and  Reynolds  [1982]  and  show  the  same  behavior.  The  behavior  of 
the  mean  flow  is  discussed  further  in  Section  5.1. 

4.5  Periodic  Velocity 

The  free-stream  excitation  imposed  on  the  test  section  and  the 
response  of  the  boundary  layer  in  all  cases  were  found  to  be  very  nearly 
sinusoidal.  Spectral  analysis  of  the  periodic  velocity  shows  that  there 
are  no  pronounced  harmonics  but  a  broad  band  distribution,  due  to  resid¬ 
ual  noise  from  the  phase-averaging  process.  The  total  harmonic  distor¬ 
tion  in  terms  of  energy  in  the  worst  case  was  about  1.5%.  Typically,  it 
was  less  than  0.5%.  Thus,  it  is  sufficient  to  present  the  amplitude  and 
the  phase  of  the  periodic  velocity  component  at  the  exciting  frequency. 

4.5.1  /tap! itude 

Figures  4.18a-b  show  the  amplitude  of  the  periodic  velocity  compo¬ 
nent  at  the  exciting  frequency  normalized  by  the  free-stream  amplitude 
at  a  number  of  frequencies.  The  abscissa  normalization  is  the  same  as 
In  the  mean  profiles.  Several  observations  can  be  made,  as  follows: 

•  At  low  frequencies,  the  oscillation  amplitude  in  the 
boundary  layer  exceeds  that  In  the  free  stream;  It  can  be 
as  large  as  twice  the  free-stream  amplitude. 

•  The  overshoot  seen  above  is  strongly  suppressed  as  the 
frequency  Increases;  the  location  of  the  peak  moves  closer 
to  the  wall. 

•  At  low  frequencies,  there  Is  a  region  in  which  the  ampli¬ 
tude  variation  Is  log-linear.  This  region  corresponds 
closely  with  the  log-linear  region  in  the  mean  velocity 


profile.  As  the  frequency  increases,  the  extent  and  the 
slope  of  this  region  decrease;  at  high  frequencies,  there 
is  no  evidence  of  such  a  region. 

•  At  high  frequencies,  the  outer-region  amplitude  is  very 
nearly  the  same  as  that  of  the  free  stream.  All  the  var¬ 
iation  is  confined  to  a  thin  region  close  to  the  wall. 

The  size  of  this  region  decreases  with  frequency.  There  is 
a  region  where  the  amplitude  is  somewhat  larger  than  the 
free-stream  amplitude. 

•  The  variation  in  the  inner  region  with  frequency  is  non¬ 
monotonic.  There  are  intermediate  frequencies  for  which, 
throughout  the  boundary  layer,  the  velocity  amplitude  is 
less  than  that  in  the  free  stream. 

Thus  one  can  conclude  that  the  periodic  velocity  field  is  very 
strongly  affected  by  the  frequency  of  oscillation.  The  behavior  at 
other  X  stations  is  similar,  as  can  be  seen  from  the  results  presented 
in  Jayaraman,  Parikh,  and  Reynolds  [19823.  Further  evidence  of  this  can 
be  seen  in  the  behavior  of  the  phase  of  the  first  harmonic  component, 
presented  next. 


4.5.2  Phase 

Figures  4.19a-b  show  the  variation  of  the  phase  angle  of  the  fun¬ 
damental  component  of  the  periodic  velocity  with  respect  to  the  phase 
angle  of  the  free-stream  velocity  at  different  frequencies.  Several 
observations  can  be  made,  as  follows: 

•  At  very  low  frequencies  (f  ~  0),  the  time  for  the  bound¬ 
ary  layer  to  respond  fully  to  the  Imposed  changes  in  the 
free  stream  will  be  much  less  than  the  oscillation  time 
period;  thus  there  will  bfe  no  phase  difference  In  the 
boundary  layer  with  respect  to  the  free  stream. 

•  At  low  frequencies,  almost  all  of  the  boundary  layer  lags 
with  respect  to  the  free  stream;  this  lag  Increases  with 
frequency  up  to  a  certain  frequency. 

•  At  low  frequencies,  rapid  changes  In  the  phase  take  place 
In  the  outer  region;  there  Is  a  region  extending  to  the 
wall  In  which  there  are  no  appreciable  phase  variations; 
this  region  corresponds  to  the  log-linear  region  of  the 
mean  velocity  and  the  amplitude  profiles  at  corresponding 
frequencies. 

•  There  Is  a  frequency  beyond  which  the  Inner  region  leads 
with  respect  to  regions  farther  from  the  wall;  the  extent 


of  the  region  of  constant  phase  in  the  Inner  region  seems 
to  decrease  and  disappear  altogether  at  high  frequencies. 

•  At  high  frequencies,  there  Is  no  phase  variation  in  the 
outer  region;  It  lags  the  free  stream  by  a  small  angle. 

Rapid  phase  variations  take  place  In  a  thin  region  close 
to  the  wall;  this  region  exhibits  phase  leads  with  respect 
to  the  free  stream.  The  phase  lead  Increases  with  fre¬ 
quency. 

These  results  for  the  phase  of  the  first  harmonic  component  confirm 
the  strong  dependence  on  frequency  of  the  periodic  velocity  field  seen 
in  the  behavior  of  the  amplitude.  Jayaraman,  Parikh,  and  Reynolds 
[1982]  give  results  at  other  X  stations  showing  the  same  trends.  The 
periodic  flow  behavior  and  its  implications  on  quantities  such  as  the 
displacement  thickness  and  the  wall  shear  stress  are  discussed  in  detail 
throughout  Chapter  5.  The  results  for  the  periodic  streamwise  turbulent 
fluctuations  will  be  presented  next. 

4.6  Periodic  Turbulence 

The  signal-to-noise  ratio  In  the  phase-averaged  turbulence 
«uV>)  Is  much  lower  than  In  the  phase- averaged  velocity.  There  could 
be  significant  harmonic  content  In  the  periodic  component  of  turbulence 
(u'u').  However,  only  the  component  at  the  exciting  frequency  is  con¬ 
sidered  In  the  following  results. 

4.6.1  Amplitude 

Figures  4.20a-b  show  the  variation  of  the  amplitude  of  oscillations 
In  <u'u'>  as  a  function  of  frequency  at  one  X  station  in  each  ex¬ 
periment.  The  amplitude  Is  normalized  by  the  square  of  the  free-stream 
mean  velocity  and  the  distance  from  the  wall  by  the  Coles  boundary  layer 
thickness  of  the  mean  velocity  profile.  Several  points  to  note  are  as 
follows: 

•  At  low  frequencies,  there  are  two  distinct  regions  where 
significant  turbulence  oscillations  take  place.  One  is  in 
the  Inner  region,  and  the  other  In  the  outer  region  of  the 
boundary  layer.  Relatively  little  oscillation  Is  present 
In  the  region  In  between. 


•  In  general,  as  the  frequency  Increases,  the  oscillation 
amplitude  decreases  throughout  the  boundary  layer. 


•  At  high  frequencies,  the  outer-region  oscillations  are 
almost  completely  suppressed,  Indicating  that  the 
turbulence  does  not  vary  through  the  cycle  and  Is 
essentially  the  same  as  the  mean. 

•  The  turbulence  structure  near  the  wall  does  vary  somewhat 
through  the  cycle,  even  at  high  frequencies.  The  varia¬ 
tion  Is  considerably  less  than  at  lower  frequencies. 

The  above  results  demonstrate  the  dynamical  nature  of  the  response  of 
turbulence  to  the  Imposed  unsteadiness.  Further  evidence  Is  provided  by 
the  phase  of  the  first  harmonic  component  given  below. 


4.6.2  Phase 

Figures  4.21a-b  show  the  variation  of  the  phase  of  the  oscillations 
In  turbulence  with  respect  to  the  free-stream  oscillations.  Several 
frequencies  In  both  experiments  are  shown.  The  following  observatlnos 
can  be  made: 

•  At  low  frequencies,  there  are  two  distinct  regions  that 
are  out  of  phase  with  each  other;  the  phase  variation 
within  each  region  is  small,  and  the  changeover  from  one 
region  to  the  other  is  abrupt. 

•  As  frequency  Is  Increased,  the  constant-phase  outer  region 
persists,  although  at  high  frequencies  the  data  are  rather 
scattered,  due  to  diminishing  oscillation  amplitudes  In 
that  region. 

•  A  thin  region  close  to  the  wall  seems  to  be  in  phase  with 
the  free-stream  oscillations  throughout  the  frequency 
range. 

•  At  medium  and  high  frequencies,  the  turbulent  oscillations 
away  from  the  wall  lag  the  free-stream  oscillations;  the 
lag  increases  with  distance  from  the  wall  and  frequency; 
there  is  no  region  of  abrupt  phase  transition,  as  in  the 
low-frequency  cases. 

The  results  for  the  amplitude  and  the  phase  of  the  fundamental  com¬ 
ponent  of  the  periodic  turbulence  clearly  show  Its  strong  dependence  on 
the  frequency  of  excitation. 


4.7  Phase-Averaged  Velocity  Profiles 

Figures  4.22a-d,  4.23a-d,  4.24a-d  show  the  phase- averaged  velocity 
profiles  every  15°  In  the  oscillation  cycle  at  X'  *  0.88  for  the  low- 
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amplitude  experiment  for  three  different  frequencies.  The  velocity  Is 
normalized  by  the  free-stream  mean  velocity  and  the  distance  from  the 
wall  by  the  mean  boundary  layer  Coles  thickness.  At  this  low  amplitude. 
It  Is  seen  that  the  profiles  are  very  similar  throughout  the  cycle. 
There  Is  a  definite  region  where  the  profiles  are  log-linear  at  all  fre¬ 
quencies. 

Figures  4.25a-d,  4.26a-d,  4.27a-d  show  the  phase-averaged  velocity 
profiles  at  the  same  frequencies  as  above,  at  X1  =  0.94  In  the  high- 
amplitude  experiment.  The  velocity  In  these  figures  is  normalized  by 
the  phase-averaged  free-stream  velocity.  Several  observations  can  be 
made,  as  follows: 

•  The  shape  of  the  velocity  profiles  undergoes  large  changes 
over  the  cycle  at  all  frequencies. 

•  At  low  frequencies,  the  boundary  layer  thickness  varies 
greatly  over  the  cycle;  the  variation  decreases  with  fre¬ 
quency.  At  high  frequencies,  the  boundary  layer  thickness 
is  essentially  constant  over  the  cycle. 

•  As  the  cycle-phase  angle  Increases,  the  wake  component  of 
the  velocity  profile  increases  while  the  near-wall  veloci¬ 
ties  decrease.  At  high  frequencies,  the  flow  near  the 
wall  Is  reversed  over  part  of  the  cycle  around  180°. 

•  When  the  near-wall  velocities  are  very  small  and/or  nega¬ 
tive,  there  Is  no  appreciable  log-linear  region  in  the 
profiles. 

•  The  size  of  the  region  of  reversed  flow  and  the  phase 
range  over  which  It  occurs  increase  with  frequency;  this 
can  be  seen  from  examination  of  phase- averaged  velocity 
profiles  at  a  number  of  high  frequencies  (not  shown). 

The  results  for  the  phase-averaged  velocity  profiles  are  in  agreement 
with  the  amplitude  and  the  phase  results  presented  earlier. 


4.8  Phase-Averaged  Turbulence  Profiles 

Figures  4.28a-h,  4.29a-h,  4.30a-h  show  the  phase-averaged  turbu¬ 
lence  profiles  every  15°  In  the  oscillation  cycle  at  X'  *  0.88  in  the 
low-amplitude  experiment  for  three  different  frequencies.  The  ordinate 
Is  normalized  by  the  square  of  the  mean  free-stream  velocity.  The  pro¬ 
files  are  very  similar.  The  near-wall  turbulence  Is  suppressed  as  cycle 
angle  Increases  from  0°. 
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Figures  4.31a-h,  4.32a-h,  4.33a-h  show  the  phase-averaged  turbu¬ 
lence  profiles  for  the  same  frequencies  at  X'  =  0.94  in  the  high- 
amplitude  experiments.  The  following  observations  can  be  made: 

•  At  low  frequencies,  the  turbulence  profile  undergoes  large 
changes  in  shape  over  the  oscillation  cycle.  The  peak 
moves  from  the  wall  to  the  outer  region  as  the  cycle  phase 
angle  increases  from  0°.  The  variation  decreases  as  fre¬ 
quency  is  increased;  at  the  highest  frequency,  the  outer 
region  is  essentially  invariant,  while  the  near-wall 
structure  varies  somewhat  over  the  cycle. 

•  At  low  and  intermediate  frequencies,  during  the  portion  of 
the  cycle  when  the  phase- averaged  velocity  near  the  wall 
is  very  small,  the  turbulence  also  is  very  small  near  the 
wall;  at  the  same  time,  the  turbulence  in  the  outer  region 
is  very  high.  Thus  the  turbulence  oscillations  in  the 
inner  and  outer  regions  are  essentially  out  of  phase,  as 
has  been  observed  before. 

Once  again,  these  results  are  in  agreement  with  the  behavior  of  the 
amplitude  and  the  phase  of  turbulence  presented  earlier. 


4.9  Computed  Quantities 

There  are  512  phase-averaged  velocity  profiles  available  over  the 
oscillation  cycle  for  each  case.  For  each  of  these  velocity  profiles, 
the  boundary  layer  99%  thickness,  the  displacement  and  momentum  thick¬ 
nesses,  and  the  shape  factor  are  calculated.  Note  that  these  are  not 
necessarily  the  same  as  the  phase  averages  of  the  same  variables.  The 

behavior  of  these  parameters  in  the  experiments  will  be  discussed  below. 

Figures  4.34a-e  show  the  variation  of  computed  thickness  and  shape 
parameters  over  the  oscillation  cycle  at  X'  ■  0.88  in  the  low- 
amplitude  experiment,  for  five  different  frequencies.  Figures  4.3Sa-e 
show  the  same  at  X'  *  0.94  for  the  high- amplitude  experiment.  The 
abscissa  is  the  difference  between  the  cycle  phase  angle  from  the  shaft 

encoder  and  the  free- stream  phase  angle  at  that  location.  The  quanti¬ 

ties  are  normalized  by  their  mean  values.  The  following  points  can  be 
noted: 

•  The  variation  of  the  computed  thicknesses  for  the  low- 
amplitude  cases  Is  essentially  sinusoidal;  the  high- 
amplitude  cases  exhibit  large  harmonic  distortions. 


•  The  variation  In  the  shape  factor  over  the  cycle  Is  small 
for  the  low-amplitude  cases;  the  high-amplitude  cases  show 
much  higher  variation,  especially  at  low  frequencies. 

•  In  general,  the  extent  of  the  variation  decreases  with 
Increasing  frequency;  at  high  frequencies,  there  Is  no 
appreciable  variation  In  the  boundary  layer  thickness. 

•  In  the  low-amplitude  cases,  the  thicknesses  are  essen¬ 
tially  In  phase  over  all  frequencies;  this  Is  definitely 
not  the  case  In  the  high-amplitude  cases. 

•  At  high  frequencies,  the  displacement  thickness  Is  nearly 
180°  out  of  phase  with  the  free-stream  velocity.  It 
should  be  noted  that,  at  very  low  frequencies,  all  thick¬ 
nesses  are  180°  out  of  phase  with  the  free  stream;  this 
phase  lag  with  respect  to  the  free-stream  Increases  as  the 
frequency  Increases.  Apparently,  for  the  displacement 
thickness,  there  is  a  frequency  at  which  Its  phase  lag 
with  respect  to  the  free  stream  is  a  maximum.  For  higher 
frequencies,  its  phase  lag  decreases,  and  at  sufficiently 
high  frequencies,  it  Is  nearly  180°. 

The  results  for  the  variation  of  the  Integral  thicknesses  over  the 
cycle  for  different  frequencies  serve  to  reinforce  further  the  strong 
dependence  of  the  periodic  velocity  field  on  frequency. 


4.10  Summary  of  Results 

The  results  presented  In  the  previous  sections  are  only  a  represen¬ 
tative  sample  from  the  very  extensive  data  set  that  has  been  created  in 
the  present  research.  Results  for  all  the  cases  studied  in  the  form  of 
plots  and  tables  are  presented  In  Jayaraman,  Parikh,  and  Reynolds 
[1982].  A  table  of  contents  of  the  microfiche  report  is  given  in  Appen¬ 
dix  A. 

The  results  shown  have  demonstrated  the  Insensitivity  of  the  mean 
flow  field  to  the  Imposed  unsteadiness,  while  stressing  the  strong 
dependence  on  frequency  of  the  periodic  flow  field.  A  rationale  for 
understanding  some  of  the  results  Is  given  In  the  next  chapter. 


Chapter  5 
DISCUSSION 


In  the  previous  chapter,  a  representative  sample  of  results  from 
the  experiments  was  presented.  The  results  show  that  the  experiments 
performed  provide  a  definitive  data  set  In  a  well-defined,  unsteady, 
turbulent  boundary  layer  over  a  wide  range  of  frequencies  and  ampli¬ 
tudes.  The  purpose  of  this  chapter  Is  to  try  to  understand  the  physical 
reasons  for  the  observed  behavior  of  the  flow.  The  emphasis  will  be  on 
major  effects  without  too  much  attention  to  very  fine  detail. 

In  summarizing  previous  related  research  on  the  subject  In  Chapter 
1,  It  was  pointed  out  that  several  unanswered  questions  remained.  Any 
answers,  suggestions,  or  clarifications  that  the  present  experiments  may 
provide  for  those  questions  will  be  looked  for  In  this  chapter.  Where 
possible,  the  results  will  be  compared  to  some  simple  analyses. 

5.1  Mean  Velocity  and  Turbulence 

The  results  presented  In  Section  4.4  clearly  show  that,  over  a  wide 
range  of  frequencies,  the  mean  velocity  and  the  mean  turbulence  are 
essentially  unaffected  by  the  type  of  free-stream  unsteadiness  Imposed 
In  these  experiments.  In  Section  2.5,  the  order  of  magnitude  analysis 
showed  that  the  governing  equation  for  the  mean  velocity  (Eq.  (2.28))  Is 
the  same  as  for  steady  flows,  as  long  as  the  free-stream  amplitude  Is 
less  than  about  30%  of  the  mean.  Thus,  Invariant  mean  velocity  would 
Imply  Invariant  u'v1  and  vice  versa.  Only  u'u'  was  measured  In  the 
experiments.  It  would  be  reasonable  to  expect  that  Its  Invariance  with 
frequency  Is  a  fair  Indication  of  Invariance  of  u'v'.  The  experimental 
results  show  that  mean  velocity  and  mean  turbulence  are  virtually  con¬ 
stant  with  frequency;  thus  they  are  self-consistent  with  and  supportive 
of  the  order  of  magnitude  analysis  that  led  to  Eq.  (2.38). 

The  low-frequency,  hlgh-ampl Itude  cases  (0.1  Hz  to  0.2  Hz)  do  show 
some  small  variation  with  frequency.  Figure  4.11b  shows  that  the 
Imposed  free-stream  mean  velocity  gradient  In  the  test  section  for  the 
0.1  Hz,  hlgh-ampl  Itude  case  Is  somewhat  lower  than  that  for  higher 
frequencies.  It  was  already  noted  In  Section  4.3  (see  Fig.  4.12b)  that 


the  rate  of  increase  of  free-stream  oscillation  amplitude  was  also  lower 
for  the  same  case.  Results  from  Section  4.9  (see  Fig.  4.35)  show  that 
the  boundary  layer  thickness  undergoes  very  large  variation  over  the 
cycle  at  low  frequencies  and  that  this  variation  is  suppressed  as  the 
frequency  increases.  Its  mean  value  is  also  larger  for  low-frequency 
cases. 


The  behavior  of  the  thickness  of  the  boundary  layer  at  low  frequen¬ 
cies  discussed  above  causes  the  imposed  free-stream  excitation  for  such 
cases  to  be  somewhat  different  from  the  higher- frequency  cases.  This  is 
believed  to  be  the  cause  of  the  small  variation  in  the  mean  velocity  and 
the  mean  turbulence  profiles  seen  in  the  low-frequency,  high-amplitude 
cases. 

5.2  Behavior  of  Design  Pressure  Gradient 

The  external  pressure  gradient  imposed  on  a  boundary  layer  exerts  a 
tremendous  Influence  on  its  behavior.  It  will  be  instructive  to  study 
the  behavior  of  the  imposed  pressure  gradient  for  the  present  experi¬ 
ments,  for  different  frequencies  and  amplitudes.  As  will  be  shown 
below,  such  an  analysis  is  very  useful  in  identifying  important  param¬ 
eters  that  govern  the  flow  and  different  flow  regimes  with  respect  to 
such  parameters.  For  the  sake  of  simplicity  and  easy  identification  of 
the  analytical  dependence  on  flow  parameters,  only  the  design  pressure 
gradient  is  dealt  with  below. 


Let 
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a  is  the  amplitude  of  oscillation  of  free-stream  velocity  at  the 
end  of  the  test  section  as  a  fraction  of  UQ. 
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p  is  the  instantaneous  static  pressure. 
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Using  Eq.  (5.1)  in  the  momentum  equation  applied  at  the  free  stream, 
together  with  the  definitions  above,  after  separating  mean  and  oscilla¬ 
tory  components,  the  following  relations  can  be  derived: 
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Ai  G  is  the  fundamental  amplitude  of  the  pressure  coefficient 
’  gradient, 

is  the  2nd  harmonic  amplitude,  and 


Is  the  phase  of  the  fundamental  component. 

The  time-averaged  free  stream  Is  clearly 

-  aX'^J  lb. 6; 

and  the  corresponding  steady,  pressure-coefficient  gradient  is  easily 
shown  to  be 

8Cps 

jjfr-  •  2a(l  -  aX* )  (5.7) 

where  ps  Is  steady  static  pressure. 

The  above  derivations  clearly  show  that  the  design  pressure  gradi¬ 
ent  Is  dependent  on  two  parameters — a  local  frequency  parameter  in  the 
form  of  a  Strouhal  number  (Sx)  based  on  the  X  distance  in  the  test 
section  and  the  local  mean  free-stream  velocity,  and  a  local  amplitude 
parameter  (aX1).  For  the  experiments  reported  here,  a  has  values  of 
0.05  and  0.25,  approximately;  Sx  ranges  from  about  0.1  to  10.0.  The 
behavior  of  each  component  of  the  unsteady  design  pressure  gradient  will 
be  discussed  below. 

The  time  mean-pressure  gradient  (Eq.  (5.2))  is  not  a  function  of 
frequency.  It  is  always  less  than  the  steady-pressure  gradient  corres¬ 
ponding  to  the  time  mean  free-stream  velocity  distribution  (Eq.  i5.7)). 
The  difference  between  the  two  Is  larger  for  higher  amplitudes. 

Figure  5.1  shows  the  variation  of  the  fundamental  amplitude  of  the 
design  pressure  gradient  normalized  by  the  mean  pressure  gradient  as  a 
function  of  frequency  parameter  Sx  and  the  amplitude  parameter  (aX‘). 
It  Is  clear  that  the  pressure-gradient  amplitude  is  essentially  a  func¬ 
tion  only  of  Sx.  The  dependence  on  the  amplitude  parameter  is  very 
weak. 

Figure  5.2  shows  the  variation  of  the  phase  of  the  first  harmonic 
pressure-gradient  oscillations  with  respect  to  the  free-stream  velocity 
oscillations  as  a  function  of  Sx  for  different  (aX* )  values.  For 
very  low  Sx  values  (l.e.,  f  «  0  Hz),  the  pressure  gradient  is  180° 
out  of  phase  with  respect  to  the  free-stream  oscillations,  since  the 
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pressure  gradient  increases  as  the  free-stream  velocity  decreases.  The 
plot  shows  that  the  pressure-gradient  oscillations  undergo  a  phase  ad¬ 
vance  as  Sx  increases.  For  very  high  Sx  values,  the  pressure 
gradient  lags  the  free  stream  by  90°.  The  dependence  on  the  amplitude 
parameter  is  once  again  much  weaker  than  the  dependence  on  Sx. 

Figure  5.3  shows  the  amplitude  of  the  second  harmonic  pressure- 
gradient  oscillations  normalized  by  the  fundamental  amplitude  as  a 
function  of  Sx  for  different  values  of  (aX ' ) .  The  harmonic  distor¬ 
tion  is  higher  for  higher  amplitudes;  it  decreases  as  Sx  increases. 
It  is  utmost  about  20%  of  the  fundamental  amplitude;  for  Sx  >  1,  it  is 
less  than  about  10%  for  all  amplitudes.  Hence,  one  can  expect  appreci¬ 
able  harmonic  distortion  in  the  variation  of  global  parameters  in  the 
high-amplitude  experiment  for  low-frequency,  downstream  X  station 
cases. 

To  summarize,  the  behavior  of  the  design  pressure  gradient  clearly 
shows  that,  for  the  present  experiments,  the  frequency  parameter,  Sx, 
is  much  more  important  than  the  amplitude  parameter.  Further,  there  are 
two  different  flow  regimes  with  respect  to  Sx,  as  discussed  below. 


Flow  Regimes 


Figure  5.1  shows  two  distinct  regions.  For  Sx  <  1,  the  normal¬ 
ized  pressure-gradient  amplitude  is  essentially  constant.  In  this 
regime,  the  effect  of  the  pressure-gradient  oscillations  will  not  be 
materially  different  from  that  at  Sx  =  U.  The  flow  can  be  expected  to 
be  governed  by  inertia.  In  general,  as  the  frequency  (and  hence  Sx) 
increases,  the  effects  of  inertia  are  to  reduce  the  amplitude  of  oscil¬ 


lations  of  the  periodic  velocity  and  to  cause  it  to  lag  increasingly 
when  compared  to  the  velocity  oscillations  at  Sx  =  U.  Ample  evidence 
for  this  behavior  is  seen  in  results  presented  in  Section  4.5. 


The  pressure-gradient  oscillations  begin  to  increase  rapidly  for 
Sx  >  1.  This  pressure-gradient  excitation  is  very  strong  compared  to 
the  mean  pressure  gradient  for  high  Sx  values.  The  increasing 
pressure-gradient  oscillations  counteract  the  effects  of  inertia.  Une 
can  expect  the  near-wall  velocities,  which  are  governed  primarily  by 
pressure  gradient,  to  vary  more  and  more  over  the  cycle  in  response  to 
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the  imposed  pressure-gradient  oscillations  as  Sx  increases.  This 
Increase  in  amplitude  of  near-wall  velocity  oscillations  should  cause 
the  wall  shear  stress  to  increase  in  amplitude.  Since  the  pressure- 
gradient  oscillations  undergo  phase  advances  as  Sx  increases,  one 
would  expect  the  near-wall  velocity  oscillations  and  hence  the  wall 
shear-stress  oscillations  to  undergo  phase  advances  as  well.  Once 
again,  results  from  Section  4.5  bear  ample  evidence  in  support  of  this. 
The  expected  wall  shear  stress  will  be  discussed  in  Section  5.6. 

5.4  Importance  of  Strouhal  Number:  Behavior  of  Displacement  Thickness 

From  the  study  of  the  behavior  of  the  unsteady  design  pressure  gra¬ 
dient,  one  can  conclude  that  Sx  is  a  very  important  parameter  for  the 
experiments  reported  here.  It  should  be  useful  to  establish  that  the 
behavior  of  global  quantities  for  all  the  cases  studied  can  be  collapsed 
into  a  single  curve  which  is  a  function  of  Sx.  It  is  doubtful  that  the 
collapse  would  be  perfect.  Definitely  there  would  be  history-  and  per¬ 
haps  amplitude- related  nonlinearities  that  cannot  be  collapsed  by  Sx. 
Still,  even  to  establish  that  Sx  is  the  dominant  parameter  for  the 
flow  would  be  very  useful,  if  only  in  reducing  the  number  of  plots  to  be 
presented. 

The  displacement  thickness  was  chosen  as  the  global  quantity  of 
interest  with  which  to  establish  the  importance  of  Sx.  From  the  phase- 
averaged  velocity  profiles,  the  variation  of  the  displacement  thickness 
was  calculated.  From  that,  the  amplitude  and  the  phase  of  the  first 
harmonic  component  were  determined  by  a  Fast  Fourier  Transform  harmonic 
analysis  procedure.  It  should  be  noted  that,  at  low  Sx  values  and  at 
most  downstream  X  stations  in  the  high-amplitude  experiment,  there  was 
considerably  higher  harmonic  content  In  the  displacement- thickness  vari¬ 
ation.  The  first  harmonic  was  always  the  most  energetic  component,  how¬ 
ever. 

Figure  5.4  shows  the  variation  of  the  amplitude  of  the  displacement 
thickness  as  a  function  of  Sx  for  all  the  cases  studied.  The  ampli¬ 
tude  has  been  normalized  by  the  mean  displacement  thickness  and  further 
normalized  by  the  relative  amplitude  of  the  free-stream  oscillations  to 
remove  any  amplitude  effects.  The  collapse  is  reasonably  good  beyond 
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Sx  *  1.  For  lower  Sx  values,  the  data  for  the  high-  and  the  low- 
amplitude  experimental  cases  seem  to  collapse  to  different  values.  This 
is  perhaps  due  to  ampl itude- related  nonlinearities  in  the  behavior  of 
the  displacement  thickness.  Note  that  the  amplitude  falls  off  rapidly 
until  it  reaches  a  value  of  about  1  around  Sx  =  3.  Beyond  that,  the 
amplitude  is  essentially  1,  as  it  should  be  for  a  high-frequency,  thin 
Stokes  layer  (see  Section  5.5). 

Figure  5.5  shows  that  the  phase  of  the  displacement  thickness  with 
respect  to  the  free  stream  collapses  reasonably  well  as  a  function  of 
Sx.  At  very  low  Sx  values,  the  displacement  thickness  oscillates  180° 
out  of  phase  with  the  free  stream.  As  Sx  increases  but  remains  less 
than  about  2,  the  effect  of  inertia  is  to  cause  the  phase  to  lag  in¬ 
creasingly  with  respect  to  Sx  *  0.  The  lag  reaches  a  maximum  of  about 
45°.  As  Sx  increases  beyond  2,  the  increasing  phase  advances  in  the 
near-wall  velocity  profile  cause  the  phase  of  the  displacement  thickness 
to  advance.  At  very  high  Sx  values,  the  displacement  thickness  is 
once  again  180°  out  of  phase  with  the  free  stream,  consistent  with  the 
Stokes  solution,  to  be  shown  later  in  Section  5.5. 

It  can  be  concluded,  from  the  above  evidence,  that  Sx  "is  a  very 
useful  and,  indeed,  the  dominant  parameter  for  the  flow  studied  here. 
The  expected  behavior  of  the  wall  shear-stress  oscillations  to  be  dis¬ 
cussed  later  lends  further  support  to  its  importance. 

5.5  High  Strouhal  Number  Behavior 

It  was  shown  in  Section  2.5  that,  for  very  high  Strouhal  numbers, 
the  governing  equation  for  the  oscillatory  component  of  the  velocity 
reduces  to  the  Stokes  equation.  The  solution  to  this  equation  is  not 
dependent  on  the  mean  velocity  field;  it  is  dependent  only  on  the  local 
free-stream  velocity  variation.  The  equation  is  easily  solved,  as  shown 
below. 


For  a  periodic  free-stream  velocity  variation  of  the  orm 


where  Ug  Is  the  time-averaged  free- stream  velocity  and  A1  y  is  the 
fundamental  amplitude  of  free-stream  oscillations.  Equation  (1.30)  can 
be  easily  solved.  The  solution  is  given  below. 
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v  is  the  kinematic  viscosity. 


In  order  to  check  the  validity  of  Stokes'  solution  for  high- 
Strouhal -number  experimental  cases  reported  here,  data  from  such  cases 
were  plotted,  along  with  the  Stokes  solution*  Cases  for  Sx  of  y.b7 
and  12.40  are  from  the  low-amplitude  experiment;  the  rest  are  from  the 
high-amplitude  experiment.  Cases  with  Sx  values  of  12.40  and  lb. 68 
are  from  X  stations  in  the  recovery  section  of  the  test  rig.  The 
first  three  or  so  points  near  the  wall  for  each  case  suffer  from  the 
near-wall  data  restriction  discussed  earlier. 

Figure  5.6  shows  the  amplitude  of  the  velocity  oscillations  in  the 
boundary  layer  normalized  by  the  free-stream  amplitude  of  oscillation. 
The  distance  from  the  wall  is  normalized  by  the  Stokes-layer  thickness 
parameter,  as  defined  above.  To  the  first  order,  the  experimental 

results  agree  quite  well  with  the  Stokes  solution. 
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Some  differences  are  worth  noting,  however.  The  recovery  section 
cases  are  somewhat  lower  than  the  cases  from  the  test  section.  It  is 
possible  that  the  boundary  layer  has  not  had  enough  length  to  adjust  to 
the  sudden  removal  of  mean  adverse  free-stream  gradient  at  the  end  of 
the  test  section.  The  peak  in  the  data  occurs  somewhat  farther  from  the 
wall  than  for  the  Stokes  solution.  Significant  changes  in  the  profile 
for  the  Stokes  solution  are  confined  to  about  five  Stokes  layer  thick¬ 
nesses;  for  the  data,  it  is  more  like  ten.  This  may  be  an  indication 
that  the  Sx  values  are  not  high  enough.  The  approach  to  the  high- 
frequency  limit  may  also  be  non-monotonic. 

Figure  5.7  shows  the  phase  of  the  velocity  oscillations  in  the 
boundary  layer  with  respect  to  the  free-stream  oscillations.  The  ab¬ 
scissa  normalization  is  the  same  as  for  the  previous  figure.  Once 
again,  to  the  first  order,  the  data  are  in  very  good  agreement  with  the 
Stokes  solution.  The  recovery  section  cases  fall  above  the  Stokes  solu¬ 
tion.  The  same  reason  cited  in  connection  with  Fig.  5.6  above  is  be¬ 
lieved  to  be  the  cause  of  this. 

Other  cases  seem  to  approach  the  solution  from  below.  In  the  outer 
region,  the  data  show  a  small  but  definite  lag  with  respect  to  the  free 
stream.  Almost  all  of  this  lag  takes  place  in  the  outer  intermittent 
region  of  the  boundary  layer  and  is  perhaps  related  to  the  dynamics  of 
that  region.  Just  as  for  the  amplitude  profiles,  the  effective  thick¬ 
ness  of  the  unsteady  region  is  somewhat  larger  for  the  data  than  for  the 
Stokes  solution. 

The  evidence  presented  above  supports  the  analysis  of  Section  2.5 
that  led  to  the  Stokes  equation.  The  conclusion  that  the  high  Strouhal 
number  results  follow  the  Stokes  solution  closely  is  in  essential  agree¬ 
ment  with  the  high-frequency,  near-wall,  unsteady  pipe-flow  study  of 
Binder  and  Kueny  [1981]. 

5.6  Expected  Wall  Shear  Stress  Behavior 

Wall  shear  stress  Is  a  very  Important  parameter  in  steady  turbulent 
boundary  layers.  It  Is  of  Interest  to  study  the  variation  of  the  oscil¬ 
lations  in  the  wall  shear  stress  in  the  experiments  reported  here.  No 
direct  measurements  of  wall  shear  stress  were  made,  it  will  be  shown 


67 


below  that  a  reasonably  good  estimate  of  the  wall  shear-stress  oscilla¬ 
tions  can  be  made  from  near-wall  velocity  data.  The  actual  values  will 
be  different  from  the  estimates.  It  Is  believed,  however,  that  Impor¬ 
tant  trends  will  be  Identified. 


Near  the  wall,  U(y)  and  u(y,t)  can  be  expanded  In  Taylor's 
series,  as  shown  below: 


where 
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The  orders  of  magnitude  of  the  terms  In  square  brackets  are  evaluated 
using  the  same  estimates  as  In  Section  2.5. 
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Clearly,  for  points  very  close  to  the  wall.  Term  3  is  always  negligible. 
Since  for  low  and  medium  Strouhal  numbers. 


6'  ~  6  , 


Term  2  is  also  negligible.  For  very  high  Strouhal  numbers, 

6'  1 
«  ~  Sl/Z 

and  hence  S'  gets  smaller  as  S  Increases.  The  contribution  due  to 
Term  2  Increases  as  S  Increases. 

It  Is  clear  from  the  above  order  of  magnitude  estimates  that  esti¬ 
mates  of  the  wall-shear  oscillations  made  from  near-wall  velocity  ampli¬ 
tude  behavior  will  be  quite  reasonable  for  low  and  moderate  Strouhal 
numbers.  For  higher  Strouhal  numbers,  the  estimates  are  expected  to  be 
lower  than  actual  values.  It  Is  also  worth  mentioning  that  the  near¬ 
wall  laser  probe  volume- averaging  problem  will  be  similar  for  both  the 
amplitude  and  the  mean  velocity.  It  Is  quite  likely  that  their  ratio  is 
more  accurate  than  the  Individual  values. 

To  see  how  reasonable  the  above  approach  may  be,  profiles  of  the 
velocity  amplitude  normalized  by  the  local  mean  velocity  were  plotted  as 
a  function  of  distance  from  the  wall.  Figure  5.8  shows  one  such  plot 
at  X'  ■  0.94  for  the  high- amplitude  experiment  at  several  frequencies. 
It  Is  seen  that,  except  at  higher  frequencies,  there  is  a  substantial 
region  near  the  wall  where  the  ratio  Is  reasonably  constant.  At  the 
higher  frequencies,  It  Is  clear  that  the  estimates  will  be  lower  than 
the  actual  values  at  the  wall. 
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In  order  to  find  an  appropriate  normalization  for  the  amplltude-to- 
mean  ratio,  the  behavior  of  the  Stokes  layer  wall  shear-stress  amplitude 
Is  Investigated  below. 

From  the  Stokes  solution,  It  can  be  easily  shown  that 


(5.14) 


(5.15) 
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Figure  5.9  shows  the  variation  of  y  as  a  function  of  Sx.  All 
cases  studied  (except  low- amplitude,  X1  =  0.39  cases  which  do  not  have 
data  close  enough  to  the  wall)  are  Included  In  the  plot.  The  Stokes 
solution  also  Is  plotted.  Several  Important  observations  can  be  made, 
as  follows: 

•  Despite  the  very  wide  range  of  frequencies  and  local  amp¬ 
litudes  and  histories  that  the  cases  represent,  the  data 
collapse  quite  well  and  show  a  very  definite  trend  with 
respect  to  Sx; 

•  For  small  Sx  values,  the  shear-stress  amplitude  de¬ 
creases  with  Increasing  Sx.  Since  there  is  no  variation 
in  the  Imposed  pressure-gradient  excitation  with  sx  in 
this  range,  the  decrease  in  amplitude  is  primarily  the 
result  of  inertia. 

•  The  shear-stress  amplitude  reaches  a  minimum  around  Sx 

of  3;  for  higher  Sx  values,  the  amplitude  actually 

Increases.  This  can  be  explained  by  the  increasing 
pressure-gradient  oscillations  that  overcome  inertia, 
especially  in  the  near-wall  region. 

•  At  high  Sx  values,  y  falls  consistently  below  the 
Stokes  solution.  There  are  probably  two  reasons  for 
this.  The  Sx  values  are  not  high  enough  for  complete 
agreement  with  the  Stokes  solution.  Evidence  for  this 
was  seen  In  Section  5.5.  Further,  It  was  noted  that,  for 
high  Sx  cases,  estimates  of  the  wall-shear  oscillation 
amplitude  from  the  near-wall  velocity  amplitude  are 
likely  to  be  lower  than  actual. 

The  phase  angle  of  the  velocity  with  respect  to  the  free  stream  of 
the  first  measured  point  was  used  as  an  estimate  of  the  phase  of  the 
wall  shear-stress  oscillations.  Several  drawbacks  of  this  estimate  need 
to  be  mentioned.  The  phase  values  at  this  point  suffer  from  the  laser 
finite-probe -volume  averaging  problem.  It  Is  believed  that,  for  the 
low-frequency  cases,  this  is  not  a  very  serious  problem,  since  the  phase 
angle  varies  very  little  In  the  Inner  region  for  those  cases.  The  high- 
frequency  cases  show  rapid  variations  In  velocity  phase  close  to  the 
wall,  and  the  estimates  seem  likely  to  be  lower  than  actual  values. 
Since  the  estimates  are  based  on  measurement  at  a  single  point,  unlike 
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the  case  of  the  Mall-shear  amplitude,  more  scatter  should  be  expected. 
Still,  the  trends  shoMn  by  these  estimates  should  be  quite  representa¬ 
tive  of  the  wall  shear-phase  behavior. 

Figure  5.10  shows  the  expected  behavior  of  the  wall  shear  stress 
with  Sx  for  the  same  cases.  Despite  all  the  reservations  expressed 
above,  the  data  follow  a  well-defined  trend  with  respect  to  Sx*  The 
following  Important  observations  can  be  made: 

•  At  very  low  Sx  values,  the  wall  shear- stress  oscllla- 
stlons  are  In  phase  with  the  free-stream  oscillations. 

As  Sx  Increases,  the  wall-shear  oscillations  Increas¬ 
ingly  lag  the  free  stream.  Once  again,  the  cause  has  to 
be  Inertia,  since  the  pressure-gradient  oscillations  are 
essentially  constant. 

•  The  shear-stress  oscillations  reach  a  maximum  lag  of 
about  17°  around  Sx  of  1.  Beyond,  the  Increasing 
pressure-gradient  oscillations  (which  undergo  phase 
advances  as  Sx  Increases)  cause  the  phase  lag  to  de¬ 
crease  and  then  cause  Increasing  phase  leads.  At  very 
high  frequencies,  the  wall  shear  will  be  45°  ahead  of  the 
free-stream  oscillations,  as  per  the  Stokes  solution. 


5.7  Log  linear  Region 

Steady  turbulent  boundary  layers  follow  the  so-called  law  of  the 
wall  In  a  region  near  the  wall,  beyond  the  viscous  layer.  If  the  phase- 
averaged  velocity  profiles  follow  the  same  law,  then  In  that  region  the 
phase  of  the  velocity  should  be  constant,  as  pointed  out  by  Simpson, 
Shlvaprasad  and  Chew  [1983].  The  appropriate  analysis  Is  as  follows: 


For  the  phase-averaged  profiles,  assume  that 
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Considering  only  the  first  harmonic  terms,  the  left-hand  side  becomes 
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The  right-hand  side  Is 
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Equations  (5.20)  and  (5.21)  Imply 
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Equation  (5.23)  Implies  that 
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The  above  analysis  shows  that  the  phase  of  the  velocity  In  the  log- 
linear  region  would  be  constant.  Further,  the  amplitude  of  the  velocity 
oscillations  would  show  a  log-linear  variation.  It  must  be  noted  that 
It  has  not  been  established  that  the  "wall  shear  velocity"  parameter  in 
the  above  analysis  Is  that  due  to  the  actual  wall  shear  stress.  Thus, 
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In  general,  there  Is  no  reason  to  expect  the  constant  phase  in  the 
region  to  equal  that  of  the  Mall  shear  stress.  If  the  frequency  is 
Included  as  another  parameter  that  affects  the  velocity  profile  in  the 
region,  the  form  of  the  velocity  profiles  should  be  as  shoMn  below. 
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Clearly,  this  form  should  reduce  to  the  steady  law  of  the  wail  for 
very  small  frequencies.  It  should  be  emphasized  again  that  there  is  no 
reason  to  expect  that  the  above  form  will  be  the  same  as  the  steady  form 
for  all  frequencies. 

The  data  shown  In  Sections  4.5  (see  Figs.  4.18  and  4.19)  clearly 
show  that,  at  low  frequencies,  there  is  a  region  of  constant  velocity 
phase.  Looking  at  the  trend  toward  the  wall  of  the  velocity  phase,  one 
might  expect  that  the  phase  of  the  wall  shear  stress  Is  also  about  the 
same  for  such  frequencies.  The  corresponding  amplitude  profiles  do  show 
a  log-linear  variation  In  the  same  region. 

As  the  frequency  increases,  however,  the  extent  of  this  region 
decreases.  Tht  wall  shear  stress  seems  to  be  ahead  of  the  log-linear 
region  In  phase.  At  high  enough  frequencies,  there  is  no  evidence  of  a 
constant  phase  or  a  log- 11  near  amplitude  region,  indicating  that  the 
steady  law  of  the  wall  Is  not  valid.  The  smooth  transition  of  amplitude 
and  phase  of  velocity  with  frequency  lends  more  credence  to  the  more 
generalized  near-wall  law  shown  above. 


5.8  Periodic  Turbulence 

The  results  presented  In  Section  4.6  (see  Figs.  4.20  and  4.21) 
clearly  demonstrate  the  strong  dependence  of  periodic  turbulence  on  the 
frequency  of  excitation.  The  results  for  the  first  harmonic  amplitude 
of  turbulence  oscillations  (Figs.  4.20)  show  that  the  estimates  of  its 
order  of  magnitude  made  In  Chapter  2  (Section  2.5.1)  are  good  over  all 
frequencies.  In  particular,  note  that  the  turbulence  amplitudes  In  the 
hlgh-amplltude  cases  (Fig.  4.20b)  are  larger  than  those  In  the  low- 
amplitude  cases  (Fig.  4.20a)  by  approximately  the  ratio  of  the  relative 
free-stream  amplitudes  in  the  two  experiments,  as  per  the  quasi-steady 
estimate  made  In  Section  2.5.1. 


The  turbulence- amplitude  profiles  show  two  distinct  regions  of 
significant  turbulence  oscillations.  The  variation  in  the  outer  region 
over  the  oscillation  cycle  is  believed  to  be  associated  with  the  varia¬ 
tion  of  the  boundary  layer  thickness  over  the  cycle.  First  of  all,  at 
low  frequencies,  the  boundary  layer  thickness  is  essentially  180°  out  of 
phase  with  the  free  stream  (see  Figs.  4.34  and  4.35),  and  so  is  the 
outer-region  turbulence  (see  Figs.  4.21).  At  low  frequencies,  both 
undergo  large  variation  over  the  cycle.  Points  very  near  the  edge  of 
the  boundary  layer  are  in  the  free  stream  for  portions  of  the  cycle  and 
inside  the  boundary  layer  at  other  times.  This  shows  up  as  turbulence 
oscillations  at  those  points.  As  the  frequency  increases,  the  variation 
in  both  boundary  layer  thickness  and  outer-region  turbulence  is  sup¬ 
pressed,  and  at  high  frequencies,  both  are  essentially  constant  over  the 
cycle  (Figs.  4.20,  4.34,  and  4.35). 

The  inner-region  turbulence  oscillations  are  also  reduced  as  the 
frequency  is  increased  (see  Figs.  4.20),  indicating  perhaps  an  inertial 
response  of  turbulence  to  the  imposed  periodic  changes.  The  observed 
Increase  in  the  phase  lag  (with  respect  to  phase  at  very  low  frequen¬ 
cies)  of  turbulence  oscillations  at  locations  away  from  the  wall  as 
frequency  increases  provides  further  evidence  for  the  inertial  response 
hypothesis.  Even  at  the  highest  frequency  studied,  however,  there  is  a 
thin  region  near  the  wall  where  there  are  significant  variations  in 
turbulence  over  the  cycle.  This  may  mean  either  that  much  higher  fre¬ 
quencies  are  needed  before  the  inner-region  turbulence  variation  will  be 
completely  suppressed  or  that  the  dynamics  of  the  response  of  turbulence 
to  imposed  excitation  cannot  be  completely  explained  by  inertia  alone. 


Chapter  6 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  following  is  a  summary  of  the  major  conclusions  of  the  present 
study  and  recommendations  for  future  work.  These  are  primarily  drawn 
from  the  results  and  discussion  presented  in  the  previous  chapters. 


6.1  Conclusions 

•  A  definitive  data  set  on  an  unsteady  turbulent  boundary 
layer  subjected  to  well-defined  oscillations  in  the  free 
stream  velocity  has  been  obtained  and  should  prove  very 
useful  in  the  evaluation  of  computational  schemes  and  as 
input  to  turbulence  models  for  unsteady  turbulent  bound¬ 
ary  layer  flows. 

•  The  time-averaged  flow  field  is  essentially  unaffected  oy 
the  imposed  unsteadiness. 

•  The  periodic  flow  is  very  strongly  dependent  on  the  fre¬ 
quency  of  excitation. 

•  The  Strouhal  number  (Sx)  based  on  the  streamwise  dis¬ 
tance  in  the  test  section  and  the  local  mean  free-stream 
velocity  Is  a  very  important  parameter  for  the  periodic 
velocity. 

•  There  are  two  distinct  regimes  of  periodic  velocity  beha¬ 
vior  with  respect  to  the  Strouhal  number,  namely,  the 
inertial  regime  (Sx  <  1)  and  the  pressure-gradient- 
dominated  regime  (Sx  >  1). 

•  The  high-frequency  periodic  velocity  results  follow  the 
Stokes  solution  very  closely. 

•  The  steady  form  of  the  law  of  the  wall  does  not  seem  to 
be  valid  for  the  phase-averaged  velocity  profiles  at  high 
frequencies. 

•  Over  the  range  of  frequencies  and  amplitudes  studied, 
there  seem  to  be  no  dramatic  interactions  between  the 
organized  oscillations  and  the  disorganized  turbulent 
fluctuations.  This  conclusion  is  specific  to  the  types 
of  data  processing  and  the  analyses  employed,  and  the 
quantities  studied  in  the  present  experiments. 


6.2  Recommendations  for  Future  Work 


The  fol lowing  are  recommendations  for  future  work  In  the 
present  experimental  facility: 

•  To  complete  the  data  set,  measurements  of  other  compo¬ 
nents  of  phase- averaged  turbulent  stresses  are  needed. 

In  particular,  measurements  of  shear  stress  (<u'v'>)  In 
the  boundary  layer  and  wall  shear  stress  (<*#>)  should 
be  made  for  a  complete  understanding  of  tne  physical 
characteristics  of  the  flow. 

•  The  experimental  apparatus  In  which  the  present  study  was 
conducted  Is  a  very  versatile  one.  Detailed  studies  on 
effects  of  different  types  of  excitation  (for  example, 
abrupt  imposition  of  an  adverse  free-stream  gradient)  and 
unsteady  turbulent  boundary-layer  heat  transfer  would  be 
very  fruitful. 

The  following  recommendations  are  addressed  to  the  unsteady  tur¬ 
bulent  boundary  layer  research  community  at  large: 

•  Unsteady,  turbulent,  boundary  layer  flow  with  separation 
Is  a  wide-open  field.  Potentially,  the  effects  of  un¬ 
steadiness  on  such  a  flow  are  much  more  Important  than 
for  attached  boundary  layers.  Careful  study  of  such 
flows  over  a  wide  range  of  relevant  parameters  (for 
example,  the  frequency  parameter)  Is  essential  for  a 
complete  understanding  of  important  phenomena,  such  as 
dynamic  stall. 

•  The  question  of  Interaction  between  the  organized  and  the 
disorganized  motions  Is  still  unanswered.  Experiments  at 
frequencies  much  higher  than  those  studied  are  perhaps 
needed  to  resolve  this  question.  Further  analysis  of  the 
present  data  from  different  perspectives  may  be  helpful 
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Flow  Legend:  TBL  *  Turbulent  Boundary  Layer, 

TPF  «  fully  developed  Turbulent  Pipe  Flow, 
TCF  *  fully  developed  Turbulent  Channel  Flow 
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a  amplitude  of  oscillation  of  the  free-stream  velocity  normalized 
by  the  mean  free  stream. 

S  Strouhal  number  based  on  the  typical  values  for  the  length 
scale  In  the  streamwlse  direction  and  the  mean  free-stream 
velocity. 
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High  Amplitude  Experimental  Cases 
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Low  Amplitude  Experimental  Cases 
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Table  4.3 

Estimated  Nth-Order  Uncertainties  at  20:1  Odds  with  Un  =  730  mm/s 


Quantity 


y/6  >  0.1 


y/6  <  0.1 
y  >  0.2  mn 


[(u'u')/U  ]  x  10* 

0 

5 

<  u  >  /u0 

0.5% 

[<  u 1 u 1  >/Uq]  X  104 

10 

*1,0^0 

0.3% 

(Vi/V  0-»5 

4° 

♦l.u  A 

lAl.u/uo:  °-25 

1° 

'u'^o^  x  1q4 

6 

.  ,  ,  lAl,„v'Uo!  60 

x  10-4 

7° 

♦l.u 'u  l  ,  2 

K.uV'V  10 

x  10"4 

32° 

1 1 

0.05  mn 

K 

5  mm 

l 

5  mn 

C  >  ,  6^ 

<  0.6  mm 

I  >  »  ®2 

<  0.5  mm 

*T 

<  2  mm/s 

Figure  3.3 


Low  Amplitude 
X' =0.94  f=l . OHz 
Y/D=0.22 


Spanwise  Uniformity 


|.000000000QO 


04b 


02b 


Low  Amplitude 
X’=0.94  f*i . 0Hz 
Y/D=0 . 22 


Figure  4.1c 


Spanwise  Uniformity 


Boundary  Layer  Thickness 


X  (cm) 

Figure  4.4a 


Displacement  Thickness 


0  0.4  0.8 

(Y/Dg) 


1.2 

Figure  4.7 


Inlet  Mean  Velocity 


00 


(cm) 


Upstream  Disturbance 

o  0.10  Hz 
+  0.50  Hz 
s  1 . 00  Hz 
x  x  2.00  Hz 


x  Xx 

X  x 


X  vX 


X  X 

xx  X  XX  X 


X  XX  V^xxx 

X 


s  S  s 

s  ss 

ss  Sgs  s  “s  %% 


S  S  s 


4.  .  +  . + 

+  ++"*"  +  +-L.+  +  ,"*■ 

+  +  ++  .  4+ 

+++++ 

o  +  n  O  OOn  ++fcKf++ 
o0o°  °°Oo°o°0  0°°  ° 

High  Amplitude  ^,QOo 

X’—  0.36 
U -733mm/s 


Y  (mm)  Figure  4.8 


V 


* 


0.8 


0.6 


Inlet  Mean  Velocity 

X  x  x  x 

y  A 
X  A 
X  A 

Xx  s  s  s 

xX  S  S 

v  x  _  s 

X  C  ®  4-  +  ■*"  + 


ry  *** 


*  *  * 
O  o  o 


0.4^g°° 


2  0.2 


High  Amplitude 
X’--0.36 
Dj.=32.  7mm 
U=733mm/s 


o  0.10  Hz 
#  0.20  Hz 
+  0.50  Hz 
s  1.00  Hz 
x  2.00  Hz 


0.4 


0.8 

(Y/Dc) 


1.2 

Figure  4.9 


Inlet  Mean  Turbulence 


High  Amplitude 
X'=-0.36 
Dc-32 .  7mm 


U=733mm/s 


o  0.10  Hz 
#  0.20  Hz 
+  0.50  Hz 
s  1 . 00  Hz 
X  2.00  Hz 


x 


Imposed  Free  Stream  Mean  Velocity 


O  *rl|cu 


lag) 


6*0 


0.5Hz 


0.8  1.2  0  0.4  0.8  1.2 

(X-XJ  /L  Figure  4.15c  (X-XJ  /L  Figure  4.15d 


Low  Amplitude  X’=0.88 


High  Amplitude 


116 


High  Amplitude  X’=0.9 4 


30 


2Q\- 


10 


Low  Amplitude  X’=0.88 

"i — r  11  1 1  hi 1 — i  i  min 1 — r  i  i  i  mi r 


x  x 


X  X 


+  0.2  Hz 
#  0.1  Hz 
o  0.5  Hz 
s  1.0  Hz 
x  2.0  Hz 


XD 

(U 

CD 


cn 

CD 

-O 


0 


-10 


3-20 


CL 

I 


□:-2o 


-20 


-20 


s 

S  S  ~  XXX 

sSsSS 


X 

s 


XxXx^xXx^ 


x  X 


*•*„ 

s 


G  0  °  0  °°G 

G  °  OG 


CP 


GQ  G  # 

Q  o  -.O  * 

°°09d  x 


*  ***  *****  *******  ****^ 


+*■ 


+  .  +  +  + 


++ 


+  + 


+  T +^+  ^+ 
+  ++-H-+r 


-20“ — -  -J — I — LU-LLLl _ I . hi  i _ i  i  i  . 


10 


-3 


10 


-a 


10 

(Y/Dc) 


-i 


10° 

Figure  4.19a 


119 


High  Amplitude 


IKE 


+  0.2  Hz 
x  0.1  Hz 
o  0.5  Hz 
s  1.0  Hz 
x  2.0  Hz 


cxxx  Xx 


Low  Amplitude 


’  *0 . 88 


High  Amplitude  X ’ =0 . 94 


124 


Low  Amp.  X ’ =0 . 88  f=0.2Hz  Low  Amp.  X’=0.88  f=0.2Hz 


Figure  4.22a  Figure  4.22b 


14K 


C\J 


CD  ^  O 

•  • 

o  o 

an/<n> 


o  o 

an/<n> 


128 


tude  X  ’  =0 . 94  f=0.2HzHigh  Amplitude 


High  Amplitude  X’=0.9 4  f=0.5HzHigh  Amplitude  X'=0.94  f=0.5Hz 


>3 


High  Amplitude  X’=0.94  f=0.5HzHigh  Amplitude  X*=0.94  f=0.5Hz 


134 


Figure  4.26c  IY/UCJ  Figure  4. 26d 


tude 


4.27d 


.88  f =0 . 2Hz  Low  Amp.  X’=0.88  f=0.2Hz 


Figure  4.28c  Figure  4. 28d 


.  2Hz 


(Y/Dg)  Figure  4.28e  (Y/Dp)  Figure  4.28f 


.88  f -0.2Hz 


(Y/DJ  Figure  4.28g  (Y/DJ  Figure  4  m 


0.88  f =0 . 5Hz 


Figure  4.29a  'T/LV  Figure  4.29b 


.88  f»0.5Hz 


.88  f»0.5Hz 


.  OHz 


(Y/D 


zhct 


Figure  4.30g  u/DcJ  Figure  4.30h 


High  Amplitude  X*-0.94  f-0.2Hz  High  Amplitude  X*»0.94  f«0.2Hz 


,oi  x  c;n/s.n>] 


,o\x  i;n/<,.n>] 


gh  Amplitude  X’-0.94  f»0.2Hz 


High  Amplitude  X,=0.94  f=0.2Hz  High  Amplitude  X*=0.94  f=*0.2Hz 


roi  x  [«n/<l.n>] 


x  i;n/s.n>] 


151 


(Y/Dc)  Figure  4.31e  (Y/Dc)  Figure  4.31f 


itude  X ’ *0 . 94  f=0.2Hz  High  Amplitude  X*«0.94  f-0.2Hz 


Figure  4.31g  vY/U,J  Figure  4.31h 


Figure  4.32b 


tude  X ’ *0 . 94  f«0.5Hz  High  Amplitude  X’«0.94  f-0.5Hz 


(Y/Dc)  Figure  4.32c  (Y/Dp)  Figure  4.32d 


*H9  ’ 


Figure  4.32e  (Y/Dp)  Figure  4.32f 


itude  X ' “0 . 94  f«0.5Hz  High  Amplitude  X’-0.94  f-0.5Hz 


tude  X ’ “0 . 94  f-2.0Hz  High  Amplitude  X’-0.94  f«2.0Hz 


Figure  4.33d 


.94  f “2 . OHz  High  Amplitude  X,-0.94  f-2.0Hz 


[<Ue>/iy,  [<Dm>/Dm],  [<D,>/D,],  [<Dj>/D2]  ,  [<H>/H] 


[ciV/ig,  [<dm>/dj.  [<d,>/d,],  kd^/dj 


62 


[<UE>/iy.  [<D„>/OJ,  [<D,>/D,].  [<Dj>/D2]  ,  [<H>/H] 


Low  Amp.  X’=0.88  f=0.5Hz 


163 


[<UE>/UE].  [<DM>/D  J  ,  [<Da>/Dj 


[<UE>/ig.  [<D„>/DJ.  [<D,>/D,].  [cO^/OJ,  [<H>/H] 


Low  Amp.  X  ’  =0 . 88  f=2.0Hz 


(deg)  Figure  4.34e 


0 


100 


200 

(0-P..J 


300  400 

(deg)  Figure  4.35a 


166 


[cU^/UJ,  [<D,^/DJ.  [<D,>/D,]  ,  [<Dt>/D,]  .  [<H>/H] 


(0~PltMf)  (deg)  Figure  4.35c 


gh  Amp.  X‘-0.9 4  f*1.0Hz 


i  i  i  i  i  i 

Free  Stream 


Delta99 


Deltal 


Delta2 


V  Shape  Factor 


100  200  300  400 

(deg)  Figure  4.35d 


(deg) 


0.01 


Displacement  Thickness  Phase 


Lay 


AD-A124  265  AN  EXPERIMENTAL  STUOV  OF  THE  DVNAHICS  OF  AN  UNSTEADV 
TURBULENT  BOUNDARV  LAVER(U)  STANFORD  UNIV  CA 
THERMO SCIENCES  DIV  R  JAVARAHAN  ET  AL.  DEC  82  TF-18 
UNCLASSIFIED  ARO-16760.  2-EG  ARO-i-88  .F/G  2074 


ea 


Figure  5.9 


Appendix  A 

CONTENTS  OF  MICROFICHE  REPORT 

Detailed  results  for  all  the  cases  studied  In  the  form  of  plots  and 
tables  are  presented  In  Jayaraman,  Parlkh,  and  Reynolds  [1982]  In  micro¬ 
fiches.  There  are  16  microfiches  In  the  report.  Their  contents  are 
summarized  below. 

A.l  Contents  of  Microfiche  1 

Microfiche  1  contains  the  boundary  layer  growth  under  zero-pressure 
gradient  conditions,  the  upstream  disturbance  during  unsteady  operation, 
and  the  Imposed  free-stream  excitation  in  the  test  section.  Table  A.1 
shows  the  data  available,  whether  the  data  are  In  the  form  of  tables  or 
plots,  and  their  locations  In  the  microfiche.  The  number  after  the 
decimal  point  In  the  locations  shown  refers  to  the  page  number  on  the 
microfiche.  The  correspondence  between  the  headings  used  In  the  tables 
and  the  present  nomenclature  1$  shown  In  Table  A.4. 

A.2  Contents  of  Microfiches  2  through  16 

Microfiches  2  through  16  contain  plots  and  tables  for  all  the  ex¬ 
perimental  cases  studied  (see  Tables  4.1  and  4.2).  Several  different 
types  of  plots  are  presented.  Table  A.2  summarizes  what  the  abscissas 
and  the  ordinates  are  In  each  plot  type  and  the  type  of  scaling  used  In 
the  abscissas.  Two  different  types  of  tables  as  described  In  Table  A.3 
are  presented.  The  nomenclature  used  In  the  plots  Is  the  same  as  the 
present  one.  The  correspondence  between  the  present  nomenclature  and 
the  headings  used  In  the  tables  Is  shown  In  Table  A.4.  For  convenience, 
the  experimental  cases  are  divided  Into  eight  groups,  as  shown  In  Table 
A.  5. 

Table  A.6  shows  the  locations  where  plots  can  be  found  In  the 
microfiches.  Table  A.7  presents  the  locations  of  the  tables.  As  shown 
In  these  tables,  there  are  two  parts  to  the  location  of  a  table  or  a 
plot.  The  part  before  the  decimal  point  refers  to  the  number  of  the 
microfiche,  and  the  part  after  the  decimal  point  shows  the  page  number 
within  the  microfiche. 
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Table  A.1 

CONTENTS  OF  MICROFICHE  1 


Table-Plot 

Zero-Pressure  Gradient  Boundary  Layer  Develop 


Free- stream  velocity:  Table 


Mean  velocity  and  turbulence:  Tables  . 


Location 


.  1.02-1.15 


Upstream  Disturbance 

•  y 

Low  amplitude,  X*  -  -1.03: 

RMS  free- stream  disturbance:  Table 

High  amplitude,  X1  »  -1.03: 

RMS  free- stream  disturbance:  Table 


High  amplitude,  X*  «  -0.36: 

Mean  velocity  and  turbulence:  Tables  .  . 


.  1.18-1.26 


Imposed  Free- St ream  Excitation  In  Test  Section 


Low  amplitude: 

Mean  and  first  harmonic  results.  Tables  . 
Mean  velocity  and  turbulence.  Plots  .  .  . 
First  harmonic  amplitude  and  phase,  Plots 
Phase-averaged  velocity.  Plots  ..... 

High  amplitude: 

Mean  and  first  harmonic  results.  Tables  . 
Mean  velocity  and  turbulence.  Plots  .  .  . 
First  harmonic  amplitude  and  phase.  Plots 
Phase-averaged  velocity.  Plots  ..... 


1.27-1.29 

1.30-1.35 

1.36-1.41 

1.42-1.53 

1.54-1.62 

1.63-1.71 

1.72-1.80 

1.81-1.98 


Table  A.2 
TYPES  OF  PLOTS 


Jm. 

Ordinate 

Abscissas 

Abscissa  Scaling 

i 

U* 

Y* 

Logarithmic 

2 

(UTUT/U^)  X  104 

(Wc) 

Logarithmic 

3 

mE) 

(*/°C) 

Linear 

4 

(ITTF/u^)  x  io4 

(*/°c) 

Linear 

5 

^,u/ai,ue) 

(t/Dc) 

Logarithmic 

6 

<P1,U  ‘  P1,UE> 

(»/Dc) 

Linear 

7 

{ai,u/ai,ue) 

t»/Bc) 

LI  near 

8 

(pl,u  -  pi,ue) 

(V/0C> 

Linear 

9 

<*i.uv/5e>  ,l°4 

(</0c) 

Logarithmic 

10 

(Pl.U'U*  -  pi,ue) 

(t/»c> 

Logarithmic 

11 

lAxo,u*/®l>  “l°4 

(»/ Be) 

Linear 

12 

(pl,0'U'  ■  pi,ue) 

(»/oc) 

Linear 

13 

[<ue>/ue],C<dw>/1I99] 

[<0^/0^,  [<d2>/d23 

[<H>/H3 

(•  *  PI,UE> 

Linear 

14 

<U>/UE 

mc) 

Logarithmic 

15 

C<u'u*>/u^3  x  io4 

(»/»c) 

Logarithmic 

Table  A.4 


NOMENCLATURE  USED  IN  TABLES 
Table  Entry  Present  Nomenclature 


A1 

Al,u 

Ale 

Ai,ue 

Ales 

Ai,ue 

AT 

Al,u'u' 

CF 

cf 

DELTA 

6 

DELTA1 

*1 

0ELTA2 

«2 

PI 

♦l,u 

Pie 

♦i,ue 

P1L 

PL 

PI' 

♦l.u'u* 

REDE LI 

w» 

REDEL2 

».v* 

Udef 

u*  -  u+ 

Umean 

u 

Ue 

u 

e 

Uo 

uo 

Uos 

<4 

U+ 

u+ 

U1  IRS 

u'u' 

Ues 

u2 

e 

UTAU 

UT 

*0 

Utr 

uriRS 

Y+ 

1+ 

185 


’  v.*  a  -  •  ’  «  '  *  -  • * 


Table  A.5 


IDENTIFICATION  OF  GROUPS  OF  EXPERIMENTAL  CASES 


Experiment  X1  Group 
Low  Amplitude  0.38  1 

0.63  2 
0.88  3 
1.13  4 


High  Amplitude 


0.39 


5 


Not*:  Location  9.46  r* t«r*  to  Psg*  46  In  Mlcroflch*  No, 


Appendix  B 

EFFECTS  OF  PROBE/WALL  RELATIVE  VIBRATIONS 


V* 


t 


Significant  errors  In  the  near-wall  unsteady  measurements  could 
result  from  periodic  variation  at  the  exciting  frequency  of  the  probe 
location  with  respect  to  the  wall.  The  following  Is  a  simple  analysis 
of  the  effects  of  such  variation  on  the  measured  near-wall  velocity. 

Let 

(y)3 

u(y,t)  =  U(y)  +  A  (y)  e  *  +  u'(y,t) 

X  ,U 


and 


y,(t) 


y  +  A, 


j(ut+*  ) 

%  n 


i.y. 


where  yp  Is  the  Instantaneous  location  of  the  probe  with  respect  to 
the  wall.  Assuming  that  the  amplitude  of  the  variation  In  the  probe 
location  Is  small  and  neglecting  higher  harmonics.  It  can  be  shown  that 


U  (y)  »  U(y  ) 
m  p 


and 


u 

m 


where 


W 


j(«ot+*  ) 

e  1,yp 


um(t)  •  u(yp»t) 


Clearly,  the  measured  mean  velocity  Is  unaffected  by  the  periodic 
variation  In  the  probe  location  with  respect  to  the  wall.  The  relative 
error  In  the  measured  periodic  velocity  Is  higher  for  locations  closer 
to  the  wall.  To  get  an  estimate  of  these  errors,  assume  that 


y 
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Further,  for  unseparated  boundary  layers,  (bU/by)  Is  a  maximum  at  the 
wall.  Thus,  for  a  worst-case  estimate, 

u  ■  - ye  *  r  4-  —  A.  e  p 


P  yp 


P  i»y. 


The  maximum  relative  error  In  amplitude  Is 

E  .  JLH 

A  A,  - 
1,t  y 
w  p 

and  It  corresponds  to  no  error  In  the  phase  angle. 
The  maximum  error  In  the  phase  angle  Is 


Ep  -  2  sin 


-1  /I 


(i  ea) 


and  It  corresponds  to  no  error  In  the  amplitude. 

For  a  given  amplitude  of  probe- location  variation,  clearly  the 
errors  are  smaller  for  adverse  pressure  gradient  boundary  layers  (small 

excitations.  In  the  present  work,  one  should  expect  the  amplitude  of 
the  variation  of  the  probe  location  with  respect  to  the  wall  to  be 
larger  In  the  high-amplitude,  high-frequency  cases,  since  the  amplitude 
of  the  oscillations  In  pressure  Is  larger  for  such  cases. 

Recall  the  LDA  probe-volume  positioning  procedures  described  In 
Section  3.3.1.  this  positioning  Is  done  under  steady  operating  condi¬ 
tions  and  Is  accurately  repeatable  to  0.025  nan.  Observation  of  the 
shadow  of  the  needle  and  the  bright  streaks  of  laser  light  on  either 
side  of  It,  during  unsteady  operation  of  the  rig,  showed  no  visible 
movement  of  the  beams  with  respect  to  the  needle.  Since  the  positioning 
repeatability  Is  about  0.025  mm,  It  can  be  concluded  that  any  movement 
of  the  beams  with  respect  to  the  needle  during  unsteady  operation  was 
less  than  0.025  mm. 

The  results  for  the  amplitude  and  the  phase  of  the  periodic  veloc¬ 
ity  presented  In  Chapter  4  show  no  anomalous  behavior  near  the  wall. 
The  low-frequency  results  show  that  the  velocity  phase  In  the  near-wall 


subjected  to  high-amplitude,  high-frequency  (large  Aj 


region  Is  nearly  a  constant,  as  it  would  be  If  the  steady  form  of  the 
law  of  the  wall  Is  applicable  to  the  phase-averaged  velocity  profiles. 
Hie  results  for  the  high-frequency  cases  are  In  essential  agreement  with 
the  expected  Stokes  solution,  as  shown  In  Chapter  5. 

Thus,  the  present  experimental  data  show  no  evidence  of  significant 
periodic  variation  In  the  LDA  measuring  volume  location  with  respect  to 
the  wall. 


